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The behavior of mesoscopic particles dissolved in a dilute solution of 
long, flexible, and nonadsorbing polymer chains is studied by field-theoretic 
methods. For spherical and cylindrical particles the solvation free energy 
for immersing a single particle in the solution is calculated explicitly. Im- 
portant features are qualitatively different for self-avoiding polymer chains 
as compared with ideal chains. The results corroborate the validity of the 
Helfrich-type curvature expansion for general particle shapes and allow for 
quantitative experimental tests. For the effective interactions between a 
small sphere and a wall, between a thin rod and a wall, and between two 
small spheres quantitative results are presented. A systematic approach 
for studying effective many-body interactions is provided. The common 
Asakura-Oosawa approximation modelling the polymer coils as hard spheres 
turns out to fail completely for small particles and still fails by about 10 % 
for large particles. 
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I. INTRODUCTION 



In colloidal suspensions the depletion interaction between mesoscopic dissolved parti- 
cles and nonadsorbing free polymer chains represents one of the basic and tunable effective 
interactions (see, e.g., Ref. 1 1| for a review). For example, adding free polymer chains 
to the solvent of a colloidal solution leads to an effective attraction between the particles 
which may lead to flocculation 0. For two individual colloidal particles or for a single 
particle near a planar wall this effective interaction can be measured even directly MM. 
In view of its importance it is surprising that for a long time the interaction between 
polymers and colloidal particles has been modelled only rather crudely by approximating 
the polymer chains by nondeformable hard spheres ||[l].|3|.[|] . 

Chain flexibility has been taken into account only more recently. Mainly the following 
two cases have been considered: (a) strongly overlapping chains (semidilute solution) 
which are described within a self-consistent field theory or within the framework of a 
phenomenological scaling theory (b) nonoverlapping chains (dilute solution) which 

to a certain extent can be modelled by random walks without self- avoidance (ideal chains) 
|TfO|-p^[. In three dimensions this latter situation is closely realized in a theta solvent ||15|| . 

Besides presenting some new results for ideal chains the main emphasis of the present 
contribution is on the generic case of a good solvent and we investigate systematically 
the consequences of the ensuing excluded volume interaction (EV interaction) |16| on 
depletion effects in a dilute and monodisperse polymer solution. The interaction of long 
flexible chains with mesoscopic particles leads to universal results which are independent 
of most microscopic details [|l5| , |r^ -19l and depend only on a few gross properties such 
as the shape of the particles. By focusing on such systems we obtain results which are 
free of nonuniversal model parameters. Due to the universality of the corresponding 
properties it is sufficient to choose a simple model for calculating these results. For 
example, in a lattice model the interaction between a particle and a nonadsorbing chain 
can be implemented as the purely geometrical restriction that the chain must not intersect 
the particle ||12|| . For our investigations we use an Edwards- type model [|l^,|17],|18[ for the 
polymer chain which allows for an expansion in terms of the EV interaction and which 
is amenable to a field-theoretical treatment. The basic elements in this expansion are 
partition functions Z[ ](r,r') for chain segments without EV interaction (as indicated 
by the subscript [0]) and with the two ends of the segment fixed at r and r'. In this 
coarse grained description the interaction of the nonadsorbing polymer with the particle 
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is implemented by the boundary condition that the segment partition function vanishes 
pi| as r or r ' approaches the surface S of the particle |H)| , |TU 



i.e., 



Z [0] (v,r') - 0, r^S. (1.1) 

The only relevant property which characterizes one of the interacting polymer chains 
is its mean square end-to-end distance 7£| in the absence of particles and other chains. 
Within the perturbative treatment of the EV interaction it will be necessary to generalize 
the three-dimensional space to a space of D spatial dimensions. In this respect it is 



convenient [17] to introduce 

K 2 X =K 2 E /D, (1.2) 

the mean square of the projection of the end-to-end distance vector onto a particular 
direction, say, the x-axis, in the D-dimensional space. For industrially produced polymers 
such as polystyrene values of TZ X up to the order of /xm are easily accessible. 

The simplest particle shapes relevant for applications are spheres and rods [l[ but the 
particles can also have more complex structures such as those of closed bilayer membranes 
in the case of vesicles |21). We note that the radius R of spherical particles can be quite 



small as compared to accessible values of 7Z X , e.g., R m 0.012 /zm in the case of Ludox 
silica particles [B3[. Rodlike objects are provided, e.g., by fibers or colloidal rods 



semiflexible polymers with a large persistence length l v such as actin for which i v « 17/zm 
24], and microtubuli [^4]. The ratio of the length I and the radius R of rodlike particles 



may be of the order of 40 or larger, in conjunction with a quite small radius such as 
R ~ 0.007 /im in the case of colloidal boehmite rods (E3L As the interaction between 
rodlike particles and polymers is concerned we consider long rods, i.e., R,7Z X <C I, and 
neglect effects which may arise due to their finite length I. In order to be able to treat 
spheres and cylinders in a unified way and in general dimensionality, we are thus led to 
consider a generalized cylinder K with an infinitely extended 'axis' of dimension 5. Such 
a generalized cylinder has been introduced in Ref. |TJ]], in the following denoted as I. The 



'axis' can be the axis of an ordinary infinitely elongated cylinder (5 = 1), or the midplane 
of a slab (6 — D — 1), or the center of a sphere (5 = 0). For general integer D and 5 the 
explicit form of K is 

K = jr = (r_L, r || ) G R D ~ S x R s ; \r±\ < i?} (1.3) 

with rj_ and rn perpendicular and parallel to the axis, respectively. Note that is a 
(i-dimensional vector with 
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d = D-5 . 



;i-4) 



The radius R of the generalized cylinder K is the radius in the cases of an ordinary 
cylinder or a sphere and it is half of the thickness in the case of a slab. For the slab 
the geometry reduces to the much studied case of (two decoupled) half spaces |19| . We 
stress that the generalization of D to values different from three is introduced only for 
technical reasons because D uc = 4 marks the upper critical dimension for the relevance 
of the EV interaction in the bulk [15,17,18|- Eventually we are interested in — and will 
obtain results for — the experimentally relevant case D — 3. These results concern the 
solvation free energy for a single particle and the depletion interaction between particles. 



A. Solvation free energy of a particle 

We consider the increase in configurational free energy of a dilute solution of long 
flexible polymers with number density n p upon immersing a single particle. For 5 > 
we actually consider a generalized cylinder with a large but finite axis length I s (i.e., an 
ordinary cylinder with axis length I or a slab with cross section area I 0-1 ) and study the 
increase n p fj>^ in free energy per fcgT and per I s in the limit I — > oo, for which I drops 



"PJ K 



out 1 25]. For a sphere n p f^ is simply the free energy increase per /c#T. The additional 
increase in free energy upon immersing the particle in the polymer free (i.e., n p = 0) 
solvent is regarded as a background term which, in an experiment, can be determined 
separately. In the asymptotic regime where both 1Z X and R are large on the microscopic 
scale (such as the monomer length or the diameter of the solvent molecules) it turns out 
that fjp takes the scaling form 

fP = R d Y d , D (x) , (1.5) 

where Y^d is a universal scaling function of the scaling variable 

x = K X /R . (1.6) 

For ideal chains (no EV interaction) and d fixed the function Y^d = yJ~ is independent 
of D (compare I where fjp was denoted as 5 fx)- Results for Y^ for d — 3 (sphere) 
and d = 2 (cylinder) have been given in Ref. [|n]] and in I. Here we calculate the scaling 
function Y^d(x) for chains with EV interaction perturbatively in terms of e = 4 — D with 
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the upper critical dimension D uc = 4. In particular we investigate the following features 

(a) For short chains, i.e., x -C 1, we assume that Y^^ix) is analytic so that it can be 
expanded into a Taylor series around x = 0. This is plausible since for short chains the 
thickness ~ 1Z X of the polymer depletion layer is much smaller than the particle radius R 
so that a small curvature expansion is applicable to fjp in which a volume term ~ R d is 
followed by a surface term ~ R d ~ l and by successive terms ~ R d ~ 2 , R d ~ s , etc., generated 
by the surface curvature. We note, however, that it can be rather difficult to actually 
prove this assumption. 

The first Taylor coefficients of the expansion of Yd^(x) around x = also determine the 
curvature energies of a particle K of more general shape provided its surface S is smooth 
and all principal radii of curvature are much larger than the polymer size TZ X (compare 
Ref. and I). Consider the increase Fjc in configurational free energy upon immersing 
a particle K, with finite volume tvc into the dilute polymer solution with bulk pressure 
n p kBT . Due to general arguments P7J in three dimensions one expects an expansion of 
the Helfrich-type 

F K - n p k B Tv K = J alS { Aa + AK X K m + Ak 2 K% + Ak g K g + . . . } (1.7a) 

s 



K m=l(^r + ^r) (l.7b) 



with the local mean curvature 

K = 

and the local Gaussian curvature 

K G = l/(RxR 2 ) , (1.7c) 

where R\ and R 2 are the two principal local radii of curvature. We use the convention 
that Ri, R 2 > means that the boundary surface is bent away from the polymer solution 



located in the exterior of K,. Provided that the expansion (|1 . 7a|) is valid the surface 
tension Act and the curvature energies Ak\, Ak 2i and Akq are determined uniquely by 
the special cases that fC is a sphere and a cylinder, respectively. Our explicit results 
for Yd.D{x) provide a strong indication that the Helfrich-type expansion ( |L7p is indeed 
valid and, moreover, does yield quantitative estimates of the surface tension and of the 
curvature energies for the polymer depletion problem in the presence of EV interaction. 
These values are the extra contributions (as indicated by the A's) to the solvation free 
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FIG. 1. Situations of short and long chains in which the limiting behavior of the scaling 
function Y ( i > d(TZ x /R) can be applied: (a) For 1Z X -C R the function Y^d determines the change 
of the surface tension Act and curvature energies Aki, Ak2, and Aug m the Helfrich-type 
expansion ( |1.7| ) of a membrane upon exposing one side of it to a dilute polymer solution, (b) 
For 1Z X 3> R the polymer can deform in order to avoid the space occupied by the particle and 
coil around a spherical (or rodlike) particle and the function Y^d exhibits the power law ( |1.8| ) 
with the Flory exponent v. 



energy of a particle in addition to its background value for the polymer free solvent (i.e., 
Up — 0), not included in Eq. ( |1.7a| ). To the best of our knowledge this is the first check of 
the expansion (|1.7|) for a nontrivial interacting system that can be realized in nature. 

For other types of systems the expansion ( |1.7|) can be violated. For example, as pointed 
out by Yaman et al. |29| , a somewhat counter-intuitive behavior arises for the case in which 
a surface is exposed on one side to a dilute solution of thin rigid rods (needles): even for 
arbitrarily small surface curvature the free energy in this case cannot be expanded in 
the analytical and local form of the Helfrich-type expansion (1.7). However, for flexible 
ideal chains instead of needles the expansion is known to apply (see I). In particular 
the asymmetry in the curvature contribution ~ R~ 2 between the in- and outside of a 
spherical or cylindrical surface as reported by the above authors for needles does not 
occur for flexible ideal chains f30 ]. 

We note that the curvature energies are experimentally accessible. For example the 
expansion (|1.7|) determines the change in surface tension and in the first and second order 
curvature energies of a flexible surface such as a membrane upon exposing one side of it 
to a solution of polymers which are depleted near the membrane (see Fig. [I] (a)). Thus 
the addition of polymers to a solution of closed membranes, i.e. vesicles, should influence 
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the phase diagram of vesicle shapes in a quantitatively controllable way (see, e.g., Ref. 
Sip. An additional experimental access to the solvation free energy will be discussed at 



the end of this subsection. 

(b) For long chains, i.e., s;>l, a single chain can deform in order to avoid the space 
occupied by the particle and coil around a spherical or rodlike particle (see Fig. 0(b)). In 
this case it turns out that Y d)D (x) exhibits a power law 

y d)D (x-oo) - A dtD xV" (1.8) 

with a dimensionless and universal amplitude A^d, provided 

d > 1/v (1.9) 

so that fjp vanishes for R — > (see Eq. ( |1.5|) ). Here v is the Flory exponent characterizing 
the power law dependence TZ X ~ N v of 1Z X on the number N of monomers per chain if 
N is large. The properties described by Eqs. ( |1.8|) and (|1.9| ) follow from a small radius 
operator expansion (SRE) (see Sec.|I B| below). 

Finally, we emphasize that fjp is experimentally accessible by monitoring the depen- 
dence of the number density n c of the colloidal particles on the number density n p of the 
polymers in a sufficiently dilute solute of immersed particles which is in thermal equilib- 
rium with a surrounding ideal gas phase with given partial pressure of the particles 



32fl . Accordingly n c is determined by a Henry- type law 

D (°) 

n c = A- 1 (1.10a) 
k B T 

where A measures the change of the solubility of the colloidal particles due to the presence 
of the polymers and is given by 

A = ex V {n p f^l s ) . (1.10b) 

For the dilute immersed particles the reduced free energy increase n p f^ I s constitutes a 
reduced one-particle potential or, equivalently, an increase in chemical potential, so that 
Eq. ( |1.10 ) follows upon equating the chemical potentials of the particles in the ideal gas 



phase and in the solution phase. 
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B. Colloidal particles with small radii 



We consider the case in which a polymer chain interacts with a spherical or cylindrical 
particle whose radius R — albeit being large on the microscopic scale — is much smaller 



than the size 1Z X of the chain and other characteristic lengths [33]. In this limiting case 
the effect of the spherical particle upon the configurations of the chain can be represented 
by a 5-function potential located at the center of the particle which repels the monomers 
of the chain. For a generalized cylinder K with a small radius R this 5-function potential 
is smeared out over its axis. Thus the Boltzmann weight Wfl-{yi} for the chain |Q arising 



from the presence of K (whose axis includes the origin) is replaced by 

W K { Yi } - 1 - A d D R d ~ x l v wk (1.11) 

with 

[ j d 5 r\\ p{r ± = 0, r ) , d < D , 
w K = < R« (1.12) 
1 P(0) , d = D , 

provided d > \ jv. The positions {y^ ; % — 1, . . . , N} of the N chain monomers which define 
the chain configuration appear in Eq. (|1.12|) in terms of the modified monomer density 



N 



i=i 



The sum of (^-functions in Eq. ( |1.13| ) is the usual monomer number density at a point r. 
We have chosen its prefactor such that p(r) is less dependent on the microscopic monomer 
structure (i.e., on what is considered as a monomer) than the sum itself. In particular 
J dPr p(r) = 7Zx^ u is independent of these details while N is not. The scaling dimension 
D — 1/v of p(r) equals its naive inverse length dimension so that the exponent of R in 
Eq. Ql.ll ) follows by comparing naive dimensions. The amplitude A^r> is dimensionless 



and universal |3 



The monomer positions {yj} are statistical variables so that Eq. fll.ll ) is a relation 



between fluctuating quantities which is to be used inside polymer conformation averages 
such as the ratio of polymer partition functions with and without the presence of K. One 
can use Eq. (|1 . 1 1|) for a variety of different situations. If K is the only particle within 
reach of the polymer chain, Eq. ( |1.11| ) leads to the free energy change given by Eq. (|1.5|) in 



the limit discussed in Eq. (|1.8|). This is the reason why the same amplitude A^d appears 
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in Eqs. ( |1.8| ) and ( p. .1 1[ ) . If there are in addition other particles or walls K', Eq. ( |1 . 1 1| ) 
can be used to calculate the polymer mediated free energy of interaction (potential of 
mean force) between K' and K (compare I and Sec. |I C| below). Equation ( |1 . 1 1|) simplifies 
the theoretical treatment of these problems significantly because K is replaced by the 
monomer density p(r). While the remaining, simpler averages depend on the particular 
problem under consideration, the universal amplitude A^d is always the same. 

In this work we study the small radius expansion (|1 . 1 1 ) for the generalized cylinder K 
for the case of polymers in a good solvent. Our main objective is to present quantitative 
estimates for the universal amplitudes A 3i3 and A 2t 3 corresponding to a sphere and to an 
infinitely elongated cylinder in three dimensions. The cylinder (i.e., d — 2) is particularly 
interesting since in this case the EV interaction changes the behavior qualitatively: while 
for ideal chains a thin cylinder is a marginal perturbation which can lead to a logarithmic 
behavior |10| and for which Eq. (|1 . 1 1|) does not apply, for chains with EV interaction the 



power law exponent d — 1/V »s 0.30 is positive and Eq. ( 1.11 ) holds. This peculiarity for 
d = 2 is reflected in the ^-expansion of A^d for D = 4 — e. 



C. Interactions between particles 

Polymer mediated interactions between particles are in general not pairwise additive, 
i.e., they cannot be written as a superposition of pair interactions []I|,|12"|. For a dilute 
polymer solution with polymer density n p we consider the total increase in reduced con- 
figurational free energy n p f t ^ upon immersing spherical particles Ki, ... , K M centered 
at r 1; . . . , r M . The quantity f^P has the form 

M M 

/t if (ri , . . . ,r M ) = £ /« + £ fS^r,) (1.14) 

pairs 
i<j 

+ ■■■ + /i£°.,K M (ri, ••• ,t m ) . 

The m-body contributions for 2 < m < M on the rhs of Eq. (|1.14j) are defined 
inductively by considering first two particles in order to define via Eq. Q1.14Q , then 
three, and so on. For spherical particles the dimension of is that of a volume, i.e., 
of (length) D . The existence of polymer mediated nonpairwise interactions has first been 
noticed within the PHS approximation, which consists in replacing the polymer by a hard 
sphere 0. Here we consider the limit for which the polymer is flexible and much longer 



than the particle radii, i.e., TZ X 3> R, and where the small-radius expansion (1.11) gives 
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a simple and quantitative description. We find that the polymer mediated interaction for 
particles with small R is drastically different from the depletion interaction for large R in 
which case the PHS approximation is reasonable and has been widely used. This confirms 
the generally accepted belief that for the applicability of the PHS approximation a large 
size ratio R/TZ X is crucial and refutes an opposite claim in Ref. [9(a)]. 

As illustration we consider three spherical particles A, B, C with radii R A , R b , Rc 
much smaller than their mutual distances and than 7Z X . It is easy to see that f tot is 
determined by the Boltzmann weights of the particles introduced in the text preceding 
Eq. fll.lip in the form 



f®(r A ,r B ,r c ) = J d D y {1 - W A W B W c } y (1.15) 

where { } y denotes the average over all conformations of a single chain in free space (i.e., 
no particles) under the constraint that one end of the chain is fixed at the point y. In 
the limit of small radii R one finds by using Eq. ( |1 . 1 1| ) that in addition to the one-body 
contributions f^\ \ and f^\ each exhibiting the scaling form described by Eq. (|1.5|) 
in the limit given by Eq. (|1.8| ), there arise two-body contributions 



f}% - - (A DjD ) 2 (R A R B ) D -^ C 2 (r A ,r B ) , (1.16a) 

f^Qi and f B 2 Q, and a three-body contribution 

fJUkc - (Md) 3 (RaRbRc) D - 1/v C 3 {r A ,r B ,r c ) . (1.16b) 

The arrows in the above relations indicate the leading behavior for small radii. Here C% 
and C3 are pair and triple correlation functions corresponding to 

C m (ri,r 2 , . . . ,r m ) = J d D y {p(n)p(r 2 ) . . . p(r m )} y (1.17) 



of the (modified) monomer density p(r) defined in Eq. ( |1.13| ) for a single polymer chain 



in free space. Since 1Z X and the relative distances r AB = \v A — *b\ are large on the 
microscopic scale these correlation functions exhibit the scaling forms 

C 2 (r A ,r B ) = TZ^- D g(z AB ) , (1.18a) 

with z AB = r AB /Tl x , and 
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C 3 (r A ,r fl ,r ) = TZ^- 2D h(z AB} z AC , z BC ) , (1.18b) 

which follow from the scaling dimension D — \jv of p(r) . Thus for three spherical particles 
with equal radii R and with center to center distances tab, tac, t bc which are of the order 
of 1Z X but much larger than R, the three-body interaction is smaller than the two-body 
interaction by a factor ~ (R/1Z x ) D ~ l l v . 

Similar fluctuation induced, not pairwise additive interactions arise between particles 
which are immersed in a near-critical fluid mixture [36|. In this case one encounters 



order parameter correlation functions instead of the present monomer density correlation 
functions. 

The small radius expressions ( |1.16| ) cease to apply — even if the equal radii R are 
much smaller than 1Z X — if some of the relative distances between the spheres become 
comparable with R. However, there are other types of short distance expansions which 
are capable to describe these latter situations. In particular we shall discuss a 'small 
dumb-bell' expansion for a pair of spheres A, B for which both R and tab are much 
smaller than the other lengths. The structure of this expansion is similar to Eq. ( |1.11| ) in 



conjunction with the lower part of Eq. ( |1.12|) , but the amplitude corresponding to Ad,d 



now depends on the ratio tab/R- We calculate this new amplitude function for the case 
of ideal chains. 

In Sec. [0] we discuss in detail the solvation free energy for a single particle. In Sec. JT| we 
consider the depletion interaction between particles. Section [TV] contains our conclusions. 
In Appendix |A| we derive the asymptotic expansions for a small and large size ratio TZ X /R 
required for Sec.||. In Appendix || we discuss the perturbative treatment of the small 
radius operator expansion. Finally, in Appendix^, we derive a short-distance amplitude 
which characterizes the behavior of monomer density correlation functions in free space 



as needed in Sec. 111. 
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II. SOLVATION FREE ENERGY OF A PARTICLE 



The free energy for immersing a particle in a dilute solution of freely floating chains 
with or without self-avoidance can be expressed in terms of the density profile of chain 
ends in the presence of the particle (compare, e.g., Eq. (3.7) in I). For the scaling function 
introduced in Eq. (|1.5|) this implies 

Y d , D (x) = ^ + n d Q d , D (ri) , 77 = x 2 /2 , (2.1) 

with tt d = 2rr d/2 /T(d/2) the surface area of the ^-dimensional unit sphere and 

00 

Qd, D (v) = J d PP d - 1 [l-M E (p,r])} . (2.2) 
1 

In Eq. ( |2.2j ) the scaling function M E (r±/ R,rj) is the bulk normalized density profile of 
chain ends at a distance r± — R from the particle surface. In Sec. [11 A| we derive the 
explicit form of Q dt D{v) i n the presence of EV interaction to lowest nontrivial order in 
e = 4 — D. In Sees. [II B| and [II C| we discuss the resulting behavior of Y dtE (x) in the limit 
of short and long chains, respectively. Finally, we obtain in Sec. |II D| approximations for 
the full scaling function Y djE (x) corresponding to a sphere and a cylinder in D = 3. 

A. Density of chain ends and polymer magnet analogy 

We employ the polymer magnet analogy (PMA) in order to calculate the density 
profile M. E of chain ends in a dilute solution of chains with EV interaction which arises 
in the presence of the nonadsorbing generalized cylinder K introduced in Eq. ( |1.3| ). As in 
I we define Ai E as bulk normalized so that it approaches one far from the particle. It is 
given by 

M E (r ± ;L ,R,u ) (2.3) 

= J d D r' Z(r,r';L ,R,u ) j J d D r' Z b (r, r'; L , u ) . 
v v 

Here Z and Zj, are partition functions of a single chain with the two ends fixed at r, r' 
in the presence and absence, respectively, of the generalized cylinder K (the subscript b 
stands for 'bulk'). The volume V available for the chain is the outer space V = M D \ K of 
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K. The parameter w characterizes the strength of the EV interaction and L determines 
the monomer content or 'length' of the chain such that 2L equals the mean square TZ% 
of the projected end-to-end distance of the chain in the absence of K and of the EV 



interaction, i.e., for w = 0. The usual arguments of the PMA |lqJ17H 19| carry over to 
the present case and imply the correspondence 

Z(t,t';L ,R,u q ) = C to ^ Lo (^(r^r ')> (2.4) 

Af=0 

between Z and the two-point correlation function ($i(r)$i(r')) in a 0(Af) symmetric 
field theory for an A^-component order parameter field <fr = ($1, ... , <J>jv) in the restricted 
volume V = M. D \ K. In Eq. ( |2.4| ) the operation 

Ao-A, = ^-J dt e Lot ° (2.5) 
c 

acting on the correlation function is an inverse Laplace transform with C a path in the 
complex t -plane to the right of all singularities of the integrand. The Laplace-conjugate 
to of L Q and the excluded volume strength uq appear, respectively, as the 'temperature' 
parameter and the prefactor of the (<fr 2 ) 2 -term in the Ginzburg-Landau Hamiltonian 

H K m = |^r{i(V$) 2 + |$ 2 + |($ 2 ) 2 } (2.6a) 
v 

which provides the statistical weight exp(— 7ii^{$}) for the field theory. The position 
vector r covers the volume V and its boundary, which is the surface of K. In order to be 
consistent with Eq. we have to impose the Dirichlet condition 

*(r) = , if |rj = R , (2.6b) 

on the boundary. This corresponds to the fixed point boundary condition of the so-called 
ordinary transition for the field theory. For our renormalization group improved 

perturbative investigations we use a dimensionally regularized continuum version of the 
field theory which we shall renormalize by minimal subtraction of poles in e = 4 — D 
(this is related via Eq. ( |2.4|) to a corresponding procedure in the Edwards model |r7|-[19|). 



The basic element of the perturbation expansion is the Gaussian two-point correlation 
function (or propagator) ($j(r) $j(r'))[ ] where the subscript [0] denotes u = 0. It is 
given by 
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($ i (r)$ J -(r')) [ o] = <%G(r,r ';*<,,#) = <% G(r ± , r[, #, |r„ - r,,'] ; * , 



(2.7a) 



5ij x < 



n=0 

oo 



(TP 



(2vr) 5 

Wt\$) G n (r x ,ri;t ,R) 



exp[t P (r,, - r„')] G n (r ± , r[; S,R) , d < D 



d = D 



\ n=0 



where a = (d — 2)/2, 5 = P 2 + to? r ± = l r _i_| ; an d $ is the angle between r_|_ and 
(compare Fig. 1 in I). Note that for d = D (last line in Eq. ( |2.7a| )) there is no parallel 
component ry — r,,' and hence no Fourier variable P. The functions Wn ($) are given by 



(2vr d / 2 )- 1 r(a) (n + a) C£(cos#) , d ^ 2 , 



(2.7b) 



(27T)" 1 (2 - 5 n>0 ) cos(ra?) 



d = 2 



where T is the gamma function, are Gegenbauer polynomials |40| , and S n fl = 1 for 
n = and zero otherwise. The functions are normalized so that J dVt d Wn = <5 n ,o • 
The propagator G n has the form 



( r (< r >h-« K. 



a+n \ 



S r 



(2.7c) 



x 



I a+n {VSR) 
K a+n (VSR) 



where r± = min(r^,r^) and r± = max(r±,r^_). For d = D the variable S is replaced 
by to. I a and K a denote modified Bessel functions [ |l| . 

The numerator in the density profile Ai e in Eq. ( |2.3| ) can be obtained from the inte- 
grated two-point correlation function, i.e., the local susceptibility x, for t > 0. Due to 
rotational invariance around and translational invariance along the axis of K the local 
susceptibility x onr y depends on the radial component r_|_ of the point r = (rj_, ry ). The 
loop expansion of x reads 



x(r±;t ,R,u ) = X m {r±;t ,R) + u oX ll] {r ± ;t , R) + O{u 2 ) 



where the zero-loop contribution is given by the integrated propagator 



(2.8) 



X m (r±;t ,R)= J d D r'G(vy;t ,R) 

v 



R 2 r 

To 



1 - 



(2.9) 
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The greek symbols on the rhs denote dimensionless variables expressed in terms of the 
radius R of K: 



r = t Q R 2 , p = r±/R 



(2.10) 



According to standard perturbation theory the one-loop contribution is given by 



UoX [1] (r±;t ,R) 

N + 2u Q , D 



(2.11) 



3 2 



N + 2 u 
3 T 



d u y G(r, y ; t , i?) G(y, y ; t , i2) x lU| (l/±; to , R) 



R 2 ^ / dtp^- 1 G(p,ip,To)g(ip,T ,e)XM(iP,To) 



where ip = y± j R (compare Eq. (|2.10|) ). The functions in the integrand of the last line in 
Eq. ( |2.11D are dimensionless and defined by 



xM^,T ) = R- 2 x [0] (y±;t ,R) 



X12) 



g(p,ip,T ) = R 2a G n=0 (r ± ,y ± ;S = t ,R) 



(2.13) 



ylg(il>,T ,e) = y[ R~ 2a G(y, y; t Q , R) 



The function g can be split into g = + g s where 



i- £/2 ^- £ r(e/2-l) 



stems from the bulk contribution of G(y, y; to , R) and 

oo 



n=0 



(2vr) 5 ' q q 



;2.14) 



(2.15a) 



(2.15b) 



K a+n (y/(f 



T 



K a+n {ip \/q 2 + T 



x < 



la+n ( / 

k K a+Tl ( ^/tq" 



[K a+n (if) y/T^ 



d< D, 



d = D 



Note that 5 = D — d = A — e — d. In the case d < D we shall consider (i (and a = (d— 2)/2) 
as a variable which is independent of .D = 4 — e whereas in the case d = D the variable 
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ft 



1 — e/2 depends of course on e. One can check that in the case d — 1, for which 
Wn C# = 0) with a = — 1/2 contributes only forn = and 1, the upper part of Eq. ( |2.15bp 
leads indeed to the half-space result 



G(y,y;t ,R) - G 6 (y,y;t ) 



(2.15c) 



>D-2 



( dP - expHVi^ + to (y±-R)], d = l, 







2V-P 2 + 1 



where the integral can be expressed in terms of a modified Bessel function. We add the 
following two remarks about the behavior of g s for d > 1 if R — > oo or R — > 0. 

(i) It is instructive to see how the behavior for the half-space arises by taking the limit 
R — > oo with to an d y± — R fixed. Consider, e.g., the case d = D = 4 corresponding to the 
sphere in four dimensions. Since upon approaching the above limit the arguments of the 
Bessel functions in the lower Eq. ( |2.15b|) become large and since many terms contribute in 
the sum over n one has to use the uniform asymptotic expansion of the Bessel functions 
(compare, e.g., sections 9.7.7 and 9.7.8 in Ref. [40(a)]) and may replace the sum by an 
integral. This yields that G(y,y;t ,R) — Gb(y,y;t ) for d = D = 4 does indeed tend 
to the half-space expression on the rhs of Eq. fl2.15cp with D = 4, where the role of the 
length P of the wavevector P is taken by the ratio n/R. 

(ii) For d > 2 and fixed nonvanishing lengths y±_ and t the quantity g s (ip, tq, e) has 
a finite limit for R — > 0, i.e., 

9 s (as) (V> V^o > £ ) = Jjm 9s OA , To , e) 



)2-<2 



(27T) 



(2.15d) 



dk k 6 - 1 (k 2 + ^ 2 r ) Q K a (y/k* + lP 2 T ) 



n 2 



7T 



d/2p, 



tt' 



k (^ 2 r ) a [K a (VV^)] 5 



d < D 



d = D 



which depends only on the i?-independent product ip^/r^ = y±\/to and describes the 

— 1/2 

behavior of g s for R ^ y±, t . This is consistent with the operator expansion for 
small radius R of the Boltzmann weight representing K when applied to a Gaussian field 
theory (compare I). While decays exponentially for ipyffo — > oo it approaches a finite 
constant for ip^7\i — > 0, which equals — a/(47r 2 ) for e = and characterizes the behavior 
of g s for R *C t This should be compared with the behavior g s ~ — (-0 — I) 2-15 

— 1/2 

which applies close to the surface of K, i.e., for < y± — R <C R, t 
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The reparametrizations |39[ 

u = 16 ir 2 f(e) pf Z u u , Z u = l + 0(u) 



(2.16a) 



and 



p 2 Z t t = p 2 ( 



2 , M + 2 u 
' 1 1 + — ^— - + 0(u 2 / 



(2.16b) 



of the bare bulk parameters Uq and to i n terms of their renormalized and dimensionless 
counterparts u and t are not affected by the presence of the surface Pl| , |38|| . Here p is 
the inverse length scale which determines the renormalization group flow and f(e) = 
1 + e fi + 0(e 2 ). The coefficient fi drops out from universal quantities and therefore can 
be chosen arbitrarily. Equation (|2.16b|) implies the renormalized counterpart 



{pR) 2 t 



(2.16c) 



of To. The renormalized, i.e., pole-free, local susceptibility Xren is related to x by [|38l ,j39H 

Xren(r±;t,R,u) = x(r±; t ,R,u ) / Z<f,(u) (2.16d) 
= X(r±] t ,R,u ) + 0(u 2 ) 



with the renormalization factor Zq> of the field <3? which deviates from one only in second 
order in u. The only pole in x^ is due to the bulk contribution gb in Eq. ( p. 15a ). When the 
results for x^ and uq x^ i n Eqs. ( |2.9| ) and ( |2.11| ) are substituted into Eq. ( |2.8D and when 
the bare parameters r and Uq are expressed in terms of their renormalized counterparts 
t and u according to Eqs. ( |2.16| ), the poles in e cancel indeed [41|. This cancellation can 
be traced back to the relation 



(2.17) 



The resulting renormalized and scaled local susceptibility X r , 
order reads 

X ren (p,r,pR,u) = X [0] (p, t) 



R 2 Xren Up to One loop 



(2.18) 



M + 2 
H — u 



^ - \n(pR) + B)r^- xW{p, r) + £ d {p, r) + 



with the nonuniversal constant 
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C E 1 , ln(47r) 
where Ce is Euler's constant, and the function 

oo 

£ d (p,r) = -8tt 2 f G(p^,r)g s ^,T,s = 0)X [0] (ij,T) . (2.20) 



Since Ed belongs to the one loop contribution and because we assume that in the last line 
of Eq. (|2.11|) the order of the -^-integration and the limit e — > can be interchanged we set 
e = in the integrand on the rhs of Eq. ( |2.20| ). This implies that in the case d = D only 



£4 enters into Eq. ( |2.18| ) (compare the remark below Eq. (|2.15b|) ). The integral on the 
rhs of Eq. Q2.20D is well-defined since the divergence of g s (ip,T,E = 0) for if) \ 1 becomes 
integrable due to the Dirichlet behavior of Q and as implied by Eq. (|2.6b|) . We also 
need the bulk value (far away from K) of the renormalized local susceptibility up to one 
loop order, which reads 

X renib (r,iiR,u) = - - ^L±l-(^-\^R) + B) + 0(u 2 ) . (2.21) 



The perturbative result ( |2.18| ) can be improved using standard renormalization group 



arguments ||38|| . Although we need only the results (|2.18|) and (p.21|) for the discussion 



of the polymer depletion problem, we note that in the asymptotic limit for which r±, R, 
and the bulk correlation length £ + for t > are large compared with microscopic lengths 
the ratio 

X ren (p, r, pR, u) I X reritb (r, pR, u) -> S^(p, 7) (2.22) 

yields a scaling form expressed in terms of the universal scaling function S^/(p, 7) with 
the scaling variables p = r±/R and 7 = _R 2 /£?. The function depends on the number 
of components of on the parameter d which characterizes the shape of K, and on 
the space dimension D. While the amplitude £ + in the bulk relation £ + = £ + t ~ V W 
is nonuniversal, the exponent v{N) is universal and depends only on M and D. The 
asymptotic scaling behavior is governed by the infrared (long-distance) stable fixed point 
for which 

u = u * = + °^ 2 ) (2-23) 

and 
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'(AO 



1 1 AT +2 ^,2, 

2 4 + 8 v ; 



(2.24) 



The bulk correlation length £ + can be defined in various ways. For definiteness we assume 
that is defined [E3] as the second moment of the two-point correlation function divided 
by 2D, which implies 

H + 2 



(£o + ) 2 = (AWr 8 ^ [l - £±f e5 + (9(. 2 )] } 



(2.25) 



with the nonuniversal constant B defined in Eq. (|2.19|) . Here the curly bracket equals 
for t — 1 and u = u* and the dependence of (£ + ) 2 011 u is contained in the dimensionless 
amplitude D t which can be expressed in terms of Wilson functions corresponding to the 
renormalization group flow of t and u [p9| , p~7| [T9|1 . When Eqs. ( [2.23| ) - (]2.25|) are combined 
with Eqs. ( |2.18 ) and ( |2.21|) one finds that X ren / X ren , b at the fixed point is indeed consistent 
with Eq. (|2.22 ) and that the scaling function Sjv - is given by 



Efr(p,7) = 7^ [01 (P,7) + ^| e 7 ^(p,7) + 0(e 2 ) . 



(2.26) 



Equation ( 2.26|) provides the general result for the bulk normalized local susceptibility of 
the magnetic analogue in the presence of K. 



The density A4 E of chain ends as defined in Eq. ( fOf ) can be related to X ren = R \ 



with Xren from Eq. ( 2.16d| ), by means of Eqs. {2A) and ( p. 16 ). The result is 
M E {r±; L , R,u ) = Z ren (p,X,pR,u) / Z reriib (X,pR,u) 



(2.27a) 



where 



Z ren (p,X,/iR,u) = C T ^ x {X ren {p,T,fiR,u)} 



M=0 



(2.27b) 



is the renormalized and scaled version of the integrated chain partition function in the 
numerator of the rhs of Eq. (|2.3| ). Here C is the operation in Eq. ( |2 .5| ) and 

(2.27c) 



A = L/{pRf = Z t L /R 2 



is the scaled counterpart of the renormalized and dimensionless chain 'length' L [17-19] 
so that Xt = Lt = L to- For large r±, L , R the end density exhibits the scaling behavior 



M E (r±;L ,R,u ) -> M E {p,r}) 



(2.27d) 
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where Me is a universal scaling function of p = r±/R and the scaling variable 



V 



2R 2 



(2.27e) 



According to our definition in Eq. (|1.2|) of 1Z E = D1Z 2 as the second moment of the 
bulk partition function Zt(r, r ') the nonuniversal prefactor r$ in the asymptotic behavior 
Til/ {2D) = r 2 L 2l/ with v = v(H = 0) has the form 



r 2 = (D L (u)r {»- 2 [l - \{B + 1 - ^) + 0(. 2 )] } 



(2.28) 



with I? from Eq. (|2.19|) . The curly bracket equals 11%/ (2D) for L = 1 and u = u* and the 
dependence of r 2 on it is contained in the amplitude Dl = 1/D t with D t from Eq. (|2.25|) 
(compare, e.g., Ref. ||19| ). Obviously rj plays a similar role as the inverse of the scaling 
variable 7 = R 2 Hi in Eq. fl2.26|) in the magnetic analogue. By using Eqs. (|2.27a|) and 
( p.27b|) and by carrying out the same steps which lead to the scaling function in 
Eq. ( p.26|) of the magnetic analogue one arrives at 



M E (p, V ) = MfW + (p, V ) + 0(e 2 ) 



(2.29a) 



where 



^2.29b) 



is the zero loop, i.e., Gaussian contribution and 



M 



[i] 



(p, V ) = C T ^{£ d (p,r)} + C T ^\^[r^-X^(p,r) + ^\} 
1 - Mf 1 (p, rj) - V j- MP (p, V )\ [ In v + C E \ + v ^ Mf (p, V ) . 



(2.29c) 



Equation ( 2.29|) provides the general result for the bulk normalized density of chain 
ends Me in a dilute polymer solution in the presence of K. According to Eq. ( |2.2| ) for 
the scaling function Y^d we only need the integrated form. The terms in Eq. ( |2.29cj ) have 
been arranged such that the p-integration in Eq. (|2.2f) can be carried out in each bracket 
separately. This leads to 



Q d , D (v) =Pf ] (v) + \Pl\rj) + 0(e 2 ) 



(2.30a) 



where 
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p "'?</) - c > J^+ii^HLX 



is the zero loop, i.e., Gaussian contribution and 

,[i], x „ , x . frlnr <9 



r— >r; 



<9r 



Lr 3 / 2 iUy7)J 



^ fo) + ^ *f (V)} [ln V + C E ]+ V ^ Pj 0] fo) 



In Eq. ( 2.30cj ) we have introduced the function 



C d {rj) = C T ^{C d {T)} 



with 



C d (r) = J dpp^Sifar) 
1 

oo 

= -8tt 2 / d^ii- 1 g s (i>,T,e = 0) [X [0] ^,r)\ 



^2.30b) 



(2.30c) 



(2.31a) 



(2.31b) 



where Eq. (|2.20 ) has been used. The functions in the integrand of the last line in 
Eq. ( |2.31b|) are given by Eqs. (|2.15b|) and Q2.12j ) in conjunction with (|2~9l) . In the case 
d = D we have to consider C±(t) only (compare the remark below Eq. ( |2.20| )). 



B. Short chains: Y d ^{x) for x — > 

The aim of this subsection is to determine the surface tension Act and the curvature 
energies Ak\, Ak2, and Akq in the expansion ( |1 . 7|) to first order in e = 4 — D by 
considering the special cases that the particle K is a generalized cylinder K with d = D, 
3, and 2. 



In accordance with the discussion in Sec. [LA] we assume that the function Y d ^(x) in 
Eq. ( |2.1| ) is analytic at x = which implies that Q d ,D {v) can be expanded into a Taylor 
series in y/rj = xj \/2. In the following we determine the first three terms of this expansion. 
The expansion is consistent with the behavior 



C d (r) = C Q r- 3/2 + Cr>r~' 2 + C 2 w r"^ + 0(r 



,(d) 



(d) -5/2 



(2.32) 
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for large r = (fiR) 2 1 of the function C^(t) in Eq. ( |2.31b| ) which we verify in Appendix |A]. 
Its form 



Qd,D (v) 



2T] 1 ' 2 r e_ 



31n2 ■ 
1 " + Co 



V 



d-l _ e 
2 4 Ll 



(2.33) 



+ 



4r/ 3 / 2 f(d-l)(d-3) 

3^ 



£ /ll 51n2 



4 V 6 



- -C 



(d) 



2 „2\ 



follows from Eqs. (|2.30| ) and (|2.31|) by inserting Eq. ( |2.32|) and the large r behavior 



r 3 / 2 K a {s/r) 



-3/2 



+ *Z± T -2 + (d-l)(d-3) __ 5/2 + G(r _ 3) 
2 8 



(2.34) 



Since Co is related to the surface tension Act it should not depend on the shape of K, i.e., 
on the value of d. Using, e.g., Eqs. (|2.31b|) , Q2.14j ), and (|2.15c| ) corresponding to a planar 
wall (i.e., d — 1) one finds 



7T 7T 

"2 + 7! 



;2.35) 



The evaluation of the coefficients C} and C 2 W in Eq. ( j2.32j) for d = D, 3, and 2 i 



(<*) 



is 



carried out in Appendix [A] by extending the method explained after Eq. ( j2.15c|) to the 
next-to-leading terms. For d = D we have to consider C±(t) only and find 



a. 



(4) 



17 157T 

~6 + ~T~ 



(2.36a) 



for d = 3 we find 



and for d = 2 



(i) 



-66 + 



(3) 



(3) 



(2) 



(2) 



8011 7T 191 v^3 



7T 



128 



17 



5 7T 

IT 



V3 



71 



551 1673 tt 40 vr 
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48 



V3 



V3 



17 5vr 

18 + T~ 



221 1791 vr 43^^ 



15 



+ 



128 



(2.36b) 

(2.37a) 
(2.37b) 



(2.38a) 
(2.38b) 
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We now determine the surface tension Act and the curvature energies Ak%, Ak 2 , and 
Akq in the expansion ( |1.7[ ). To this end we need to generalize this expansion to be 
applicable to (D — l)-dimensional surfaces of general shape with values of D different 
from three. According to differential geometry for integer D > 3 the expansion has again 
the form ( 1.7a ) and the corresponding curvatures are given by [13 



i=l 



and 



K, 



G 



D-1 

£ 



RjRj 

pairs J 



(d-l)(d-2) 1 
2 & 



(2.39b) 



where Ri are the D — 1 principal local radii of curvature. The last expressions on the rhs 

of Eq. ( |2.39|) apply to the surface of a generalized cylinder K with integer d < D. These 

expressions hold because the surface of K has d — 1 finite local radii of curvature Ri = R 

which allow for (d — l)(d — 2)/2 different pairings. Note that for D = 3 Eq. ( |2.39| ) reduce 

to Eqs. (|1.7b[) and ( |1.7c|) . Applying Eqs. ( |1.7aD and (|2.39| ) to generalized cylinders K in 

D dimensions one infers from the definition (|1.5| ) of Yd,o and Eq. ( p.l|) the general form 

d-1 1 
~R 



n p k B T RQ d ,D{v) = Act + Aki 



(2.40) 



Ak 2 



{d 



+ o(fi- 3 ) 



of Q for small 77. Explicit results for Act, Aki, Ak 2 , and Akq follow from the results 
(p|) - (pi) for the coefficients C} d) by comparing Eq. ( |2.40|) with Eq. (p.33|) . Using 
7] = TZ^/(2R 2 ) we find for the surface tension to first order in e = 4 — D 



Act = n p k B TTZ x J- jl - | 



1 



31n2 



+ C 



} + 0(e' 



(2-41) 



n p k B TTZ x 0.798 (1 - 0.0508 e) + 0{e z 



Here and in the rest of this subsection by taking e = 1 one obtains the corresponding 
estimate for the physical dimension D — 3. By setting d — 2, 3, and .D in Eq. ( |2.40|) , in 
which the generalization of d to noninteger values is obvious, we find for the curvature 
energies 



Akx 



n p k B T^f{l- £ -C} 2) ) +0{e 2 ) 
n p k B TUl 0.5 (1-0. 131 e) + C(e 2 



(2.42) 
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AK 2 = -n^ B T^{l-|[(H_^)_ 8C f)]} +<D (e>) (2.43) 
« -n p k B TKl 0.133 (1- 0.0713 e) + C(e 2 ) , 



Ak g = -Ak 2 - £n p fc B T— ^=-^- + 0(^ 2 ) (2.44) 

3 v 27T 2 



and finally 

« n p k B TKl 0.133(1 -0.177 e) + C(e 2 ) . 

Note that Aacj is fixed by considering only one of the cases d — 2, 3, and .D (we chose 
d = 2 in Eq. (ggg) ). However, since A«i must not depend on the value of d one derives 



the two conditions 



/7 (3) p (4) 

Ci ( 2) = = (2 45a) 



which must be fulfilled if the expansion ( |1 . 7|) is consistent up to one loop order in the EV 
interaction of the polymer chains. Similarly, Ak 2 and Akq are fixed by considering only 
two of the cases d = 2, 3, and D (we chose d = 2 and 3 in Eqs. ( |2.43| ) and ( |2.44| )). Thus 
one derives the third condition 



C 2 



(i) 



r (3) _ r (2) 
u 2 u 2 



(2.45b) 



By using the values of C} and C 2 as derived in the cases (a), (b), (c) above one finds 
that all three conditions fl2.45| ) are indeed fulfilled. This confirms to first order in e the 
assumption preceding Eq. (|2.32| ) that the scaling function Y^d (x) is analytic at x = 



and that the Helfrich-type expansion ( |1.7|) is applicable to the present polymer depletion 
problem for chains with EV interaction. Considering the involved analytical means which 
were necessary to derive the coefficients (see Appendix |A]) we regard this as a very 
valuable and important check of our calculation and in addition as a strong evidence that 
the above statements for Y^d (x) are general properties in D = 3 which hold beyond the 
present perturbative treatment. 

Note that the EV interaction of the polymer chains reduces the absolute values of 
the surface tension Act and of the curvature energies Aki, Ak 2 , and Akq as compared 
to ideal chains. This trend can be anticipated because the EV interaction of the chain 
monomers effectively reduces the depletion effect of the particle surface (compare, e.g., 
Ref. fim). However, the corresponding corrections are relatively small so that the overall 
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behavior is changed only quantitatively. Thus exposing one side of a flexible membrane 
to a solution of polymers which are depleted by the membrane favors a bending of the 
membrane surface towards the solution and leads to a weakening of its surface rigidity. 
The sign of the Gaussian curvature energy Akq will generally favor surfaces with higher 
genuses (see the Introduction and I). If the resolution of an experimental setup is high 
enough to observe these effects quantitatively, the corrections due to the presence of the 
EV interaction of the polymer chains as compared to the behavior for ideal chains should 
be detectable. Specificly we consider the experiments for vesicles reported by Dobereiner 



et al. [^T|]. The intrinsic spontaneous curvature energy K\ of the bilayer membrane is to 
be identified with their quantity —2kc /Ra (compare Eq. (9) in Ref. [j3l] ). The difference 
A^i (see Eq. ( |2.42| )) should be added in the presence of polymers in the solution. The 



length Ra is of the order of the size of the vesicle. Upon inserting the values k ~ 10 19 J 
and Co ~ 10 (compare Fig. 9 in Ref. [31]) one infers K\Ra ~ —2 x 10~ 18 J. On the other 



hand, for T = 300 K and n p TZ^ of order unity, which means that the polymer solution is 
still in the dilute regime so that the result ( |2.42| ) is valid, one has AkiTZ x k2x 10~ 21 J. 



The size ratio TZ x /Ra is of the order of 1/100 <C 1 for realistic values Ra ~ 10 /im and 
1Z X ~ 0.1 /im. We conclude that Aki can reach a value up to about 10% of Ki in a 
quantitatively controllable way. This can be expected to lead to observable effects near a 
shape transition of the vesicle. 

Ideal chains lead to the behavior that all contributions in curly brackets on the rhs of 
the expansion (|1.7a|) of second and higher order in the curvature vanish for the case of a 
generalized cylinder with d = 3 and D > 3 arbitrary (compare, e.g., Eqs. (3.9) and (3.11) 
in I). This encompasses, in particular, the three-dimensional sphere for which d = D = 3. 
For the contribution of second order in the curvature the reason is a combination of the 
general property Kq = for d = 3 (compare Eq. (|2.39|)) with the property Akq = — Ak 2 



valid for any dimension D if the chains are ideal. However, the last property is rather 
special and is violated for polymers with EV interaction in D slightly below 4 since 
Eq. Q2.44D implies 

t-> 3 n (3) 

Ak 2 + Ak g = -en p k B T—^^— + 0{e 2 ) (2.46) 

3 v Itx 2 

w -en p k B TKl 0.0141 + 0(e 2 ) . 

There is no reason to believe that this violation is removed in D = 3. Rather the crossover 
to a behavior Q 3j3 ~ {1Z X / R) 1 ^ 11 for 1Z X /R — > oo with the Flory exponent v ~ 0.588 (see 
Eqs. (|1.8j ), ( |2.1| ), and Sec. |11 C| ) implies infinitely many nonvanishing terms in the small 



25 



curvature expansion ( |1.7a| ) in D = 3. Thus in the physically important case of the 
three-dimensional sphere the appearance of the EV interaction does lead to a qualitative 
change. 

As an illustration, consider a spherical membrane in the dilute polymer solution with 
both sides of the membrane exposed to the polymers. In this case the contributions to 
AKiK m in the expansion ( 1.7a| ) from each side cancel and Eq. ( [2.46 ) implies that for chains 
with EV interaction the free energy cost for immersing the spherical membrane is smaller 
as compared to a flat membrane with the same area. This is different from the behavior 
for ideal chains for which the solvation free energies for a spherical and a flat membrane 
with same area are equal in this case. 



C. Long chains: Y^d{x) for x —* oo 

Figure |2] shows in the (d, _D)-plane the dashed line d = l/u(D) It separates 

generalized cylinders K which are relevant perturbations for long polymer chains with 
EV interaction (such as the strip in D = 2 or the plate in D = 3) from those which 
are irrelevant and for which Eq. Ql.llQ applies. The latter are located in the shaded 
region above the line and comprise the disc in D = 2 and the sphere and the cylinder 
in D = 3 and are of main interest here. For the sphere and the cylinder in D = 3 we 
show within an expansion in e = 4 — D that the first order result for Y^d{^) given by 
Eqs. fl2.30| ) and ( j2.1|) is consistent with the expected power law (|1.8| ) and we determine 



the corresponding universal amplitude A^d to first order in e for d = D, 3, and 2. 
These results in conjunction with the known value for A 2t2 in D = 2 are used in order 
to derive improved estimates for A S3 and A 2j3 corresponding to a sphere and a cylinder, 
respectively, in D = 3. 

The line d — 1/ v(D) itself corresponds to marginal perturbations leading to a behavior 
which in general is different [47| from Eq. (1.8). We shall discuss neither this nor the 



crossover from marginal to power law behavior which may arise for points close above the 
line. Instead, in the case d = 2 and D < 4 we shall obtain the e-expansion of A 2 d by 
analytic continuation in d from the corresponding value for d > 2. 

In the following we set d to an arbitrary value with 2 < d < D. By inserting Qd,n(v) 
from Eq. (fcgOD in Eq. Q one finds 



Y d n(x — > oo) — ► VL a 



2ar, - £ -C d (r,) + 0(e 2 ) , // = .r~/2 . (2.17) 
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2 3 4 D 



FIG. 2. Diagram of generalized cylinders K which behave — in the renormalization group 
sense — as relevant or irrelevant perturbations for nonadsorbing polymers. The parameter 
d < D characterizes the shape of K and D is the spatial dimension (see Eq. Ql,3| )), The point 
(d, D) = (2, 2) corresponds to a disc in D = 2 and the points (3, 3) and (2, 3) to a sphere and an 
infinitely elongated cylinder in D = 3, respectively. The line with D = D uc = 4 and arbitrary 
d represents the upper critical dimension where the polymers behave like ideal chains and from 
which the perturbative expansion in e = 4 — D starts in order to study the effects of the EV 
interaction. The open circles indicate points (d,D) for which d — l/i>(D) = 0. These points 
are connected by the dashed line so that within the shaded region above it the power law ( |1.8| ) 
applies and K represents an irrelevant perturbation. The paths indicated by the arrows are 
discussed in the main text. 



where a = (d — 2)/2 > 0. The first term in square brackets stems from Pj (r/) in 
Eq. (|2T30bD and C d {r]) is given by Eq. (|273l|) . Both the term Q d /d on the rhs of Eq. (|2~lD 
and the sum of the terms following —Cd{rj) on the rhs of Eq. ( |2.30c|) are subdominant to 
the leading behavior in Eq. ( |2.47| ). According to Appendix [A| this leads to 

In rj s 



Y d ,D(x -> oo) -> Vt d 2ar] 
= Vl d 2a 
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(2.48) 



27 



The constant is given by 



E A 



4tt 5 



a 



ft, 



3 
2 



In 2 



*(d/2) 



(2.49) 



where for d = D we have to consider i? 4 only. The corresponding numbers ftj are: 

1 



B A = , ft, « 0.01047 , Bo 



8tt 2 



(2.50) 



The result in the second line of Eq. 
with the power law ( |1.8| ) since 



W) for the behavior of Y^pix — > oo) is consistent 



1/(21/) = 2 -l/(2^) x l/^ 



I £ In 77 
~ 4 ~2~ 



+ 0(5 2 



(2.51) 



(see Eq. ( |2.24| ) for M = 0). The universal amplitude A^d is determined by Eqs. ( |2.48| ) 
and (|2.49| ) to first order in e = 4 — D with the results 



.4 



D.D 
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1 - 21nvr 
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E 



(2.52a) 



19.739 (1 - 0.625 e) + 0(e 2 



A 



3,D 



2ttU + 



1 ln2 Ce 

2 2 2 



8ir 2 B 3 + ^ + ^ + 



+ 0{e 2 ) 



(2.52b) 



6.283 (1 + 0.490 e) + 0(e 2 ) , 



A 2tD = e2n 3 B 2 + 0{e 2 ) w 0.785 e + 0(e 2 



(2.52c) 



where Eq. ( fj.5U| ) has been used. 

From Eq. ( |2.52c ) it is evident that A 2 ,d vanishes in the limit D /* 4 which reflects 
the fact that for ideal chains, for which 1/u — 2 and the condition (|1.9|) is violated, the 
power law ( p..8|) does not apply |^| . However, we succeeded in calculating the amplitude 
A 2) d for D < 4 to first order in e = 4 — D by following a path in the (d, -D)-plane 
which circumvents the point (2,4) as indicated by arrows in Fig.||] and along which the 
power law (|1.8|) does apply with a positive amplitude A^d- Accordingly, first one has to 
exponentiate Eq. ( |2.48| ) with respect to e for a > fixed in order to obtain the power law 
( |1.8| ), and then one has to perform the limit d — 2 = 2a \ for the resulting amplitude 
A d D for D = 4 — e fixed. 
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We note that the values for A 33 and A 2t3 which follow from Eqs. ( |2.52| ) by setting 
e = 1 are estimates which depend on the path taken. For e = 1, e.g., Eqs. Q2.52ap and 



( p.52b|) lead to the different estimates 7.39 and 9.36, respectively, for the same quantity 
A 33 (the corresponding paths in the (d, _D)-plane are indicated by the two upper arrows 
in Fig. 0). This discrepancy is caused by the present perturbative calculation of A djD . 

This unpleasant feature can be cured. As mentioned in Sec.[T^| the power law 



is a special consequence of the small radius expansion (SRE) in Eq. ( 1.11 ). Via the 
polymer magnet analogy this operator expansion is related to a corresponding SRE in a 
field theory. This allows one to understand not only the mechanism behind the SRE in 
terms of perturbative field theoretic methods for D slightly below 4 (as demonstrated in 
Appendix [FJ) but also to use nonperturbative methods for D = 2 [E?| which incorporate 
the result A 2t2 = 3.81 (see the end of Appendix ||). Improved estimates for the amplitudes 
^3 and A 2>3 can be deduced by combining the e-expansion of A^d in Eq. (|2.52|) with the 
above value for A 2y2 . To this end we assume that A^d is a smooth function of d and D. We 
consider the following interpolation schemes [^UJ for the functions f(e = 4 — D) = Ap^, 



A 2 ,d, Ad-i,,D) and Aq-d,d, which appear as curves in the A^d -surface shown in Fig.|3]: 
(a) pure polynomial: 

f(e) = f tt (e) = /(0) + me + a 2 e 2 , (2.53a) 



(b) (1,1) - Pade form: 



f{e) = f b {e) = /(0) + -^f- . (2.53b) 

1 + b 2 e 



For Ad,d and A 2) d the coefficients on the rhs of Eqs. (|2.53 ) are fixed by Eqs. ( |2.52a|) 



and ( |2.52c|) , respectively, in conjunction with /(2) = A 2j2 . Note that the corresponding 
paths in the (d, D)-plane are straight lines, i.e., in particular smooth paths, so that A^d 
behaves smoothly as function of e along these paths (see Fig.|^). We obtain estimates for 
A 33 and A 2t3 by the corresponding mean values f m (e) = (f a (e) + /&( £ ))/2 for e — 1 and 
use the difference between the two values / a (l) and /&(1) as an estimate for the error. For 
the sphere this leads to 



and for the cylinder to 



A 3<3 = 9.82 ±0.3 (2.54) 



A 2t3 = 1.23 ±0.2 . (2.55) 
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FIG. 3. The universal amplitude A^d corresponds to a two dimensional surface over the 
base plane (d,D) (compare Fig||). The full dot corresponding to ^2,2 and the thick solid lines 
represent the known parts of this surface. The solid parts of the dashed and of the dotted 
lines indicate the slopes of these lines at their end points D = 4 according to Eq. ( (2.52|) . The 
dashed lines themselves including the desired estimates for ^3 and ^3 (open squares) display 
the corresponding mean values f m {e) = (f a (s) + /b(e))/2 of the two interpolation schemes 
described in Eq. ( |2.53 ). The same holds for the dotted lines and for the two values of A^d for 
(d, D) = (2.5, 3.5), which have been calculated for a self-consistency check. These two values are 
connected by the short full line in order to indicate the deviation caused by the fact that the 
two dotted lines miss each other slightly (for the exact surface A^^d, of course, the two dotted 
lines do intersect at this point). The smallness of the deviation underscores the reliability of the 
interpolation scheme. 
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So far the ^-expansion of in Eq. ( |2.52b| ) has not been used. Now it can serve as 
a check for the reliability of the interpolation method leading to Eqs. fl2.54| ) and ( 2.55 ). 
In combination with the known curve A dA = 2n d / 2 /r((d - 2)/2) Eq. (jpgg ) determines 
the plane tangent to the A dD -surface at (d,D) = (3,4) which leads together with the 
value for A 2 ^ in Eq. ( |2.55| ) to approximations of the form (|2.53|) for the curve A D _ 1D . 
Corresponding approximations for the curve A e _£, ^ follow from the known tangent plane 
at (d,D) = (2,4) and the value for A^ 3 in Eq. ( |2.54| ). The resulting mean values f m (e) 
are shown as dotted lines in Fig.|| A satisfactory self-consistency check for the accuracy 
is provided by the observation that at the particular point (d,D) = (2.5,3.5) at which 
the two exact dotted lines should cross the approximate ones in Fig. ^ are only slightly 
off by the small amount of 0.3. 



D. The complete scaling function Y d ^(x) 

The full scaling function Y d D (^) describes the crossover between its analytic behavior 
for x = TZ X /R — > and the power law ( |1.8|) for x — > oo which have been discussed in 
Sees. |lllj| and |11 Q respectively. Here we present estimates for the complete functions 
^3,3 i x ) an d Y 2: 3 (x) corresponding to a sphere and a cylinder, respectively, in D = 3. The 
global behavior of Y d ^(x) is conveniently characterized in terms of the function 



1 



d 



QdAv) 



Qd,o(x) = - Y dtD (x) - -i = n d , n = x 2 /2 , (2.56) 



x 



where Qd,D(v) is defined in Eq. ( |2.2|) . According to Eq. ( p.40|) the value 0d,i)(O) is related 



to the surface tension Act in the Helfrich-type expansion ( |1.7| ) and the first and second 
derivatives of Q dt D(x) at x = are related to the corresponding first and second order 
curvature contributions, respectively (compare Sec. |llB| ). In the opposite limit x — > oo 
the function Q dy D {%) exhibits the power law 

6 d , D (x^oo) -> A d , D x x ' v - x (2.57) 



as implied by Eqs. ( |2.56| ) and (|1.8| 



We derive estimates for 0^,3 (x) for d = 3 and 2 by a combination of (a) an appropriate 
exponentiation of the one loop order result for ®cL,d{x) for large x with (b) an interpolation 
to the regular behavior of @d,D{x) f° r small x. The exponentiation is necessary because 
the actual forms of Eqs. ( |2.30| ) and ( |2.56|) do not exhibit the power law ( 2.57|) (compare 
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Sec. [11 Cp . By construction our estimates incorporate the improved estimates for A 33 
and A 2t 3 as given by Eqs. ( 12.54] ) and ( |2.55p in conjunction with the best available value 
v(D — 3) ~ 0.588 [39(b)] for the power law (|2.57|) . For x — > they reproduce the regular 
behavior as implied by Eqs. (|2.40| ) - (|2.44|) with e = 4 — D set to one. 

(a) Exponentiation: There are numerous possibilities to add on the rhs of Eq. ( |2.30a|) 
higher order terms in e such that the power law ( |2.57| ) is reproduced. For the sphere we 



choose to consider the path d = D (compare the derivation of Eq. ( 2.54 )) and define 



[1 + 0.179 e 2 )Q D 



W( V ,e) 



x 



exp 



1 

2 - - 

v 



(2.58) 



with P 4 and p} from Eq. (|2.30|) . The superscript oo refers to the behavior for x — > oo. 
Here 

1 [^o(v^)] 2 ' 



2r2 [^(v^)] 2 . 



;2.59) 



is the expansion of Pf, (r/) to first order in e. The rhs of Eq. (|2.58|) is consistent with 
Eq. ( |2.30|) in conjunction with Eq. (|2.56| ) and leads to the power law ( |2.57|) . The constant 
0.179 e 2 has been introduced in order to incorporate for e = 1 the improved estimate 
^3,3 ~ 9.82 given by Eq. (|2.54| ). For the cylinder a proper exponentiation is more involved 
than for the sphere (compare Sec. |11 (J| ). In this case we consider 



X 



(1 + 0.412e) - 77 exp 



2-- )E 2+£ 
v 



(2.60) 



exp 



The constant Ed is given by Eq. ( |2.49|) . The rhs of Eq. (|2.60|) is consistent with Eq. (|2.30|) 
in conjunction with Eq. ( ]2.56| ) and leads to the power law ( p.57| ). The constant 0.412 e 
has been introduced in order to incorporate the improved estimate A 2 ^ ~ 1.23 given by 
Eq. ( p.55|) . In order to carry out the exponentiation procedure described by Eq. (|2.48|) we 
use the properties [|I§ sK a+1 (s) = 2aK a (s) + sK a _i(s) and K £ = K + O(e 2 ) of the mod- 
ified Bessel functions K a and K a+ i which appear on the rhs of Eq. ( |2.30b| ) (compare the 



discussion after Eq. ( ^.52c| )). For e = 1 the quantities E 2+£ , P 2+£ , and K £ / 2 in Eq. ( |2.60| ) 
correspond to d — 3. This should be compared with Eq. ( ]2.58|) for the sphere in which P4 
corresponds to d — 4. Thus in a certain sense the exponentiation procedure implied by 
Eq. ( ]2.60|) partially circumvents the point (d, D) = (2, 4) in the (d, D)-plane by connecting 
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(3, 4) with (2, 3) on a straight line in the same way as Eq. ( |2.58| ) connects (4, 4) with (3, 3) 
(see Fig.0). 

(b) Interpolation: Since OcL,d{x —* 0) behaves regularly no exponentiation is needed in 
this limit. We are thus led to introduce a function 0^°^,(x) by the first two contributions on 
the rhs of Eq. ( |2.30a| ) in conjunction with Eqs. (|2.56 ). For consistency with the discussion 



in Sec. |11 ij| (compare Eq. ( |2.4(J| )) for the sphere and the cylinder we consider the function 
Qd,D for fixed d = 3 and d = 2, respectively. So far for both the sphere and the cylinder 
we have constructed two functions: the one with superscript oo describes well the limit 
x — > oo whereas the one with superscript describes well the limit x — * 0. We smoothly 
interpolate between these two functions using the switch function 

tanl) ( C|.r — 

x 



six I = - 
v ' 2 



tanh (v!X - ^) + 1 ( 2 - 61 ) 



with vi,V2 > so that s(x — > 0) — > and s(x — > oo) — > 1. By using this switch function 
we constitute the estimate 

e d , 3 (z) = [1 - *(*)] Q^(x) + s(x) B$(x) (2.62) 

for d = 3 and 2 in D = 3. The parameters i>i,i>2 in Eq. ( |2.61| ) should be adjusted such 
that 0,2,3 (x) behaves as smoothly as possible in the whole range of x. It turns out that 
within the corresponding region of Vi,v% the relative changes of Qd,3(%) are quite small. 
We shall use the values v± = 0.1 and v<i = 1.5. Note that the rhs of Eq. ( |2.62j ) reflects 
both the behavior Q^(x — > 0) and the behavior Q^\x — > oo) since s(x) exhibits an 
essential singularity in both limits. 

The resulting functions &d,3{x) for d = 3 and 2 are shown in Fig.|]. The functions 
which enter &jfl(x) and Q d °^\x) in Eq. fl2.62|) have been derived in Sec. [II A| and we have 
carried out the inverse Laplace transforms in Eqs. ( |2.30| ) and ( |2.31|) numerically (see also 



Table [TI| in Appendix [A]) . Figure |] shows the behavior both for chains with EV interaction 
and for ideal chains. 



It is evident that the power law (|2.57|) does not only determine the asymptotic behavior 



of the scaling function for x — > oo but it also influences the behavior down to values of 
x of order unity. This implies that for a quantitative analysis it is indispensable to take 
the behavior ( |2.57| ) into account, in particular accurate values of the amplitudes A 3i3 and 
A2,3- For the cylinder and chains with EV interaction the approach towards the power 
law (|2.57 ) is rather slow which is consistent with the fact that in this case the exponent 



d — l/vtt 0.30 in Eq. ( |1.11| ) is positive but small. Note that the functions Qd,3( x ) exhibit 
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FIG. 4. Scaling function 0^3 (x) (solid lines) for (a) d = 3 and (b) d = 2 corresponding 
to a sphere and a cylinder in D = 3, respectively (see Eqs. ( |1.5[ ) and ( |2,56| ), where O3 = 4-7T 
and O2 = 27r). The lines labelled 'EV correspond to chains with EV interaction and the lines 
labelled 'ideal' to ideal chains. The dotted lines display the quadratic polynomial in x which 
characterizes the behavior Od,3(% — ► 0) (see Eqs. ( 2.40 ) and ( 2.56[ )). For a sphere and ideal chains 
©3,3(3^) is simply a linear function of x (compare Eqs. (3.9) and (3.11) in I). The dashed lines 
display the power law ( pT57|) - For the curves corresponding to chains with EV interaction the 
amplitudes Ad$ from Eqs. ( |2.54| ) and ( |2.55| ) and the value v = 0.588 have been incorporated. For 
a cylinder exposed to ideal chains the scaling function ©2,3(2; — ► 00) diverges as x/lnx instead 
of a pure power law. 



smaller values for chains with EV interaction than for ideal chains. This is consistent 
with the exponent l/i/-l« 0.70 for chains with EV interaction being smaller than the 
exponent 1/u — 1 = 1 for ideal chains. This difference in behavior is in accordance with 
the general observation that the EV interaction effectively reduces the depletion effect of 
the immersed particle (compare the related discussion at the end of Sec. |II B| ). 



III. DEPLETION INTERACTION BETWEEN PARTICLES 



First, we consider the effective interaction between a thin rod and a planar wall con- 
fining the polymer solution. This is another example which demonstrates the importance 
of the qualitative difference between the behavior for ideal chains and chains with EV 
interaction which we have discussed in Sec. [TB| . Then, we consider the effective inter- 
action between two or three small spherical particles in the unbounded solution. When 
R is small compared with 1Z X and the distances between the particles, the small radius 
expansion ( |1 . 1 1| ) applies. On the other hand if both R and some of these distances are 
small compared to 1Z X and the remaining distances, operator expansions slightly more 
complicated than Eq. ( |1 . 1 1| ) are expected to hold. In particular we shall consider a 'small 
dumb-bell' expansion for two spheres. Finally, we compare our results with those of the 
PHS model §. 



35 



A. Interaction of a thin rod with a planar wall 



In view of the depletion driven adsorption of colloidal rods onto a hard wall p0| it 
is of interest to consider a cylinder with radius R and length I immersed parallel to and 
at a distance T> of closest approach surface-to-surface from a planar wall W in a dilute 
polymer solution (compare I). We consider the special case R <C T>,1Z X and T>,1Z X <C /. 
Using Eq. (4.19) in I we obtain the corresponding effective free energy of interaction in 
three dimensions, 

AF depl (V) = -n p k B TA 2;i lR 2 (TZ x /R) 1/u [1 - M { J\v/TZ x )] , (3.1) 



with the number density n p of the polymers in the bulk solution and the bulk normalized 
density profile M^\z/lZ x ) of chain monomers in the half-space (without the cylinder) 
as function of the distance z from the wall W . This universal density profile can be 
determined experimentally, e.g., by neutron reflectivity [^T|] (compare also Fig. 5 in I). 
Note that Eq. ( [3.1|) , in which the universal amplitude A 2 ^ enters (see Eq. ( j2.55| )), is only 
valid for chains with EV interaction. Equation Q3.1|) gives rise to an attractive interaction 
between the rod and the wall. The rhs of Eq. ( |3.1| ) is fixed by well-defined quantities and 
is independent of nonuniversal model parameters |52 . 




FIG. 5. Angular and positional coordinates # and T> for a rod near a wall. The magnitudes 
R -C TZ X -C I of length scales for which Eq. (|3.4| ) is expected to hold is shown schematically. As 
in Ref. [50] we assume that the rod has two caps with radius R and that its cylindrical part has 
a length I. 
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Consider now a solution of rods whose number density is so low that the interaction 
between the rods is negligible and the rods behave like an ideal gas. In this case the 
density of the rods near the wall is determined by the one-particle potential energy 

U(V,0) = U w (V,ti) + U depl (V,$) , (3.2) 

where Uw equals infinity for rod configurations which overlap with the wall W and equals 
zero otherwise and Ud ep i is the polymer-induced effective free energy of interaction. Both 
contributions depend on T> = h — R with h the distance between the center of the rod and 
the wall and on the angle •& between the rod and the wall (compare Ref. |50| and Fig. |5|). 



The density of rods c(T>) with arbitrary orientations compatible with the presence of the 
wall is given by 

1 

c(p) = c b ^J d(costf) exp[-U(p,ti)/k B T\ , (3.3) 

-i 

where c& is the number density of the rods in the bulk solution. Even for R -C TZ X <C I, 
in general one cannot identify Ud ep i on the rhs of Eq. (|3.2|) with AFd ep i because Eq. ( |3.1|) 
holds only for a long, thin rod with R -C Tt X) T> <C I and which is in addition parallel 
to the wall, i.e., $ = tt/2. However, since the attractive effective interaction Ud ep i(T> , i?) 
vanishes beyond a scale given by 1Z X one expects that c(T>) will be largest in a region 
T> < 1Z X so that an approximation which is valid in this region should be accurate in 
general |[50|| . Due to 1Z X <C I the rods must be located almost parallel to the wall in this 
region so that it should be a good approximation to replace U ( i ep i(T>, , &) by its value for 
•d = tt/2, i.e., by AFd ep i(V) in Eq. ( j37[l) . In this case the integration over $ in Eq. (|3.3|) 
can be carried out explicitly |50| so that 



c(V) w c (V) ex P [-AF depl (V)/k B T] (3.4) 

with cq{V) = 2cbV/l for < V < 1/2 and cq(T>) = Cb otherwise. Equation (|3.4|) is 
expected to be valid provided R 7Z X ,T> <^ I. Up to a certain degree of accuracy 
it is even possible to drop the restrictions on T>: since the prefactor Cq(T>) on the rhs 
of Eq. ( |3.4| ) vanishes linearly with D, the density of rods c(T>) will be much smaller for 
T> w R than for the intermediate region T> w 1Z X where Eq. ( |3.4| ) applies. On the other 
hand, for V > I the density of rods c(V) will be close to c {V) due to 1Z X /, which is 
also consistent with Eq. ( |3.4| ) . Therefore we expect that Eq. (|3.4|) is applicable even in the 
whole range of V provided R <C TZ X <ti I. 
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B. Depletion interaction between spherical particles 



In Eqs. ( |1.16| ) - ( |1.18| ) the interaction between small spherical particles is expressed in 
terms of the universal small sphere amplitude Ad,d and the monomer density correlation 
functions C m of a polymer chain in unbounded infinite space. Numerical values of the 
former for several spatial dimensions D are summarized in Table |. For the latter we note 
the relations 



J d D r A C 2 (r A ,v B ) = TZ^ 



(3.5) 



and 



J d D r c C 3 (v A ,r B ,r c ) = TZ^C 2 (v A ,r B ) 

10 n 



(3.6) 



which follow from the defining Eqs. (|1.13|) and ( 1.17 ). Simple limiting behaviors arise if 
the relative distance r AB = \r A — r B \ — albeit being large on the microscopic scale — is 



much smaller than other mesoscopic lengths. For the pair correlation |p~5| , pT7| , p~8| , |53[1 this 
limiting behavior takes the form 



C 2 (r A ,r B ) ar AB 



tab < T^a 



(3.7) 



For the triple correlation one finds 

C 3 (r A ,r B ,r c ) -> or 



-(D-l/u) 



AB 



Co 



*a + r B 



(3. 



tab <7^,|r c -(r A + r B )/2| . 



Here a is a universal bulk amplitude. For ideal chains a = is only defined for spatial 
dimensions D > 2 for which 



a 



(id) 



71 



-D/2 T (D 



V 2 



1 



(3.9) 



For chains with EV interaction, however, a remains finite down to D = 1. Numerical 
values of a for several D are summarized in Table |. In Appendix we show how these 
values can be obtained. 
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TABLE I. Numerical values of the small sphere amplitude ^4_d,d and of the short-distance 
amplitude a for chains with EV interaction and for ideal chains. 



D 


4 


3 




2 


1 


A D ,D 

A (id) 
D 


19.739 
19.739 


9.82 ±0.3 
6.283 





3.810 
(marginal) 


(marginal) 


a 


0.101 
0.101 


0.13 
0.318 


oo 


0.278 

(marginal) 


1 



For ideal chains the correlation functions G-i and C3 can be calculated in closed form. 
For the pair correlation one finds for arbitrary D 



Ci' ld \r A ,r B ) 



(3.10) 



x 



D 

2 



with T the incomplete gamma-function and g 2 = rJ B /(27l 2 ) = z^ B /2. For D = 3 
Eq. ( |3.10|) reduces to 



C} td) (r A , r fl ) = vr^ 2 (K x /z AB ) S(z 2 B /2) , D = 3 



-3/2 



(3.11) 



where 



S(g 2 ) = (1 + 2g 2 ) v ^ei{cg - 2gexp(-g 2 



(3.12) 



is the Fourier transform of the Debye scattering function JTB] JTT| , [T5[ ] . For the triple corre- 
lation one finds in D = 3 

S((z BA + z Ac f/2) 



C3 d) {r A ,r B ,r c ) 



2vr 5 / 2 



z BA Z A c 



(3.13) 



S{{z AB + z BC f/2) S((z AC + zcb) 2 /2) 



z AB z BC 



Z A C Z CB 



D = 3. 



One can verify that the expressions ( p.llp and ( |3.13p obey the short distance relations 
(p|) and (U) with a {ld) = tt" 1 from Eq. (gg). 

The limiting behavior ( |1.16| ) ceases to apply if the mutual distance between the small 
spheres becomes comparable with the order of their radii. As an illustration we consider 
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two spheres A and B with equal radii Ra = Rb — R- While for R <C tab^x the 
reduced free energy of interaction n p f} 2 ^ is given by Eq. ( |1.16a| ) and, in particular, for 
R < r AB < TZ X by 

fi% - " (A D>D f a i^- 1 ^) nl'» r A ^ lM (3-14) 
due to Eq. (|3.7|), one finds for R, tab *C with arbitrary tab/R that 

/i? - - (2 - M) Ac, D R D - 1 '" Tl x J v . (3.15) 

Here 

M = M(V/R) (3.16) 

is independent of TZ X and is a universal function ofV/R with 

P = r AB -2i? (3.17) 

the distance of closest approach surface-to-surface between the two spheres A and B. 
Equation ( |3.15| ) holds because on the large length scale set by 1Z X the 'dumb-bell' com- 



posed of the two spheres with small R and T> can be considered in leading order |j4] as 
a pointlike perturbation as in Eq. ( |1.11| ) in conjunction with the lower part of Eq. ( 1.12 ), 



but with the amplitude A^d replaced by an amplitude- function M.Ad,d which depends 
oiaV/R [j55] . Consistency of Eqs. ( |3.15|) and (|3.14j ) requires that 

2 - M -»• A DtD a {V/R)-^- 1 ^ , V/R -> oo . (3.18) 

In the opposite limit T>/R — > of two touching spheres [19(b)] the function A4(T>/R) 
approaches a constant larger than one because the dumb-bell is a stronger perturbation 
than a single sphere. Similar to the small radius expansion ( 1.11 ) the small dumb-bell 



expansion has — via the polymer magnet analogy — its counterpart in the A/"-component 
field theory. An easy way to obtain the explicit expression for A4Ad,d is to calculate the 
energy density profile (— 3> 2 (r))db,crit at the critical point of the field theory in the presence 
of the two spheres with radius R and Dirichlet boundary conditions which represent the 
dumb-bell (db) centered at the origin, and to compare the result with the corresponding 
result as derived from the small dumb-bell expansion in the form 

A A A T„ ( — ^ ( r ))db,crit ~~ ( — $ )fe crit / Q 1 n \ 

MA D n = hm — nn 1 / / t / r\\ ^o/ w • (3.19) 

r-oo RD-V" (#(0) $ 2 (r)) 6iCrit 
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Here \I/ is the normalized energy density introduced in Eq. ( |U7|) and the rhs of Eq. ( |3.19| ) 
is taken in the limit Af \0. The evaluation of the numerator is simplified by means of a 
conformal transformation relating it to the corresponding quantity between two concentric 
spheres [36(b)]. 

For ideal chains — corresponding to a Gaussian field theory — the latter quantity is 
known and leads to 

M = 2 - e 1 ' 2 ) - 2 jr ( D ~ 3 + l ) [Q-^ D - 2 )/ 2 ) + , (3.20a) 

i=o V 1 ' 

where 9 is related to the dumb-bell parameter T>/R via 

1 . « i \ „V 1 /£> x 2 



Equation ( |3.20|) provides an explicit expression for fj^ in Eq. ( p.!5|) . In particular one 



can check Eq. (|3.18| ) by using the relations A^%a^ = 2 and v 1 = 2 valid for ideal 



chains. In D = 3 one finds the leading behavior [56 



(2) 



I A 



which determines not only the solvation free energy of but also the depletion force between 
two small touching spheres in a dilute solution of ideal chains in D = 3. Numerical eval- 
uation of Eq. ( |3.20|) for arbitrary V/R in D = 3 shows that the crossover of /ab/^/^) 
from the behavior given on the rhs of Eq. ( |3.21| ), valid for T> R TZ X , to the behavior 



—R/T>, valid for R -C T> -C TZ X , is monotonic and without inflection point. Since this 
holds also for the crossover from 7Z X to R 1Z X CDas implied by Eqs. (|3.11 



and ( |1.16a|) , one finds that upon increasing the distance T> the reduced free energy of in- 



(2) 

teraction n p / A ^ between two small spheres is monotonically increasing and the attractive 
force ^npf^B is monotonically decreasing in the whole range of T>. 

This is different from the behavior of a particle with small radius interacting with a 



planar wall (compare Sec. [Ill A| and I). In this case the attractive force -^AF dep i is not 



monotonically decreasing with increasing T> but exhibits a maximum at a distance T> max of 
the order 1Z X since the monomer density profile in Eq. ( |3.1| ) has a point of inflection. 

This qualitatively different feature applies not only for a thin cylinder but also for a small 
spherical particle near a wall ||57|| . Another remarkable difference between the two cases 
is the behavior of the force in the limit R, T> 1Z X . While in the case of two spheres 
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~§D n pfA^B increases as TZ^ V for 1Z X — > oo, the force -^AF dep i between the particle and 
the wall exhibits a finite limit for 1Z X — > oo. This is plausible since the particle eventually 
moves into a region which is already depleted due to the presence of the wall. 

It is interesting that for two touching spheres in a solution of ideal chains in D = 3 
the form of the normalized interaction free energy 

fi%/(2R) 3 = - a - {TZ x /Rf , R<t:n x , (3.22a) 



for small radius as implied by Eqs. (|3.21|) , (|L5 ), ( |1.8|) and A3 = 2tc is very similar to 
its counterpart 

fi%/(2R) 3 = - h - (U x /R) 2 , R^>n x , (3.22b) 
for large radius following from the Derjaguin approximation | 5q| , which is supposed to 



be exact in this limit. Both forms display the same power in the length ratio TZ X /R and 
their amplitudes a = tt(1 — In 2) = 0.964 and b = (tt/2) In 2 = 1.09 are nearly the same. 
Although we do not have an explicit expression for the normalized interaction free energy 
for 1Z X /R of order unity we expect that either of the two limiting forms (|3.22j ) provides 



a reasonable approximation even in the intermediate regime. This is confirmed by the 



computer simulation results of Ref. |12| in which the chain is modelled as an iV-step 
random walk on a simple cubic lattice with N = 10 or 100 and the diameter 2R of each 
of the two touching spheres equals 10.5 lattice constants. This corresponds to the values 
0.06 or 0.60 of \{1Z X /R) 2 = |A^/(|l0.5) 2 and each of our two forms leads to estimates 
which are fairly close |59| to the simulation results 0.04 or 0.50 displayed in Fig. 3 of Ref. 



13 for the quantity -VL 2b */a* ol ^here which is to be identified with —f^ B /{2R) here 



In order to be able to appreciate the results for the depletion interaction of particles 
with small radii R as obtained in this subsection it is instructive to compare them with 
those of the PHS model || extrapolated to the case of small R | 60| . A force displaying a 



maximum at a distance T> max of order TZ X for the effective interaction between a particle 
and a planar wall and a monotonical decrease of the force with increasing T> for two 
particles of equal size are also found within the PHS model when extrapolated to R <C 1Z X . 
However, the PHS model does not produce the decrease of the absolute values of the free 
energy of interaction and of the force with decreasing R but in the limit R — ► rather leads 
to finite quantities which are independent of R. For example, in the case of a thin cylinder 
or a small sphere near a wall the maximum force in the PHS model is not proportional 
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to R d x l v IZx^ 1 as for a flexible chain but rather to 1Z£ 1 . In the particle- wall case the 
PHS model also fails to predict that the force becomes independent of 1Z X for R, V -C 1Z X . 

Even for the much studied case of a large sphere radius, i.e., R 3> TZ X , for which the 
PHS approximation is expected to work best, the deviation of the PHS approximation 
from the Derjaguin result is considerable. The PHS approximation implies b = 1 in 
Eq. ( |3.22b|) i.e., it leads to a free energy f\ B for two large touching spheres of equal 
size whose absolute value is too small by about 10%. The same ratio (vr/2) In 2 between 
the Derjaguin result and the PHS approximation appears in the case of a single large 
sphere touching a planar wall (compare footnote [28] in Ref. [|62| ). 



IV. SUMMARY AND CONCLUDING REMARKS 

We have studied the interaction of mesoscopic particles (spheres, cylinders, and planar 
walls) with a dilute solution of long, flexible, free, and nonadsorbing polymer chains which 
are depleted by the particles in good or theta solvents. The properties for a single particle 
as well as the effective interaction between two or more particles have been considered. 

One topic of main concern has been to investigate in a systematic and quantitative 
way how the excluded volume (EV) interaction between the chain monomers modifies the 
ideal chain behavior. Our results are in line with the plausible conjecture that weaker 
depletion effects arise from chains with EV interaction than from ideal chains with the 
same Flory radius. Another main topic has been the description of situations in which 
the particle radius R is small compared with the Flory radius 1Z X so that the chain will 
coil around the particle (compare Fig.|l|) and in which the classic PHS treatment ignoring 
chain flexibility || is clearly of no use. For example, consider the limit R/TZ X — > in which 
the spherical or cylindrical particle degenerates to a point or a thin needle, respectively, 
on the scale of 1Z X : for flexible polymers both the solvation free energy of the particle and 
its polymer mediated free energy of interaction with other particles vanishes in this limit 
whereas these two quantities remain finite for the rigid polymers of the PHS model. 

Our analysis is based on the polymer magnet analogy which maps the polymer problem 
with interactions within a single polymer chain and between a polymer chain and a particle 
onto a Ginzburg-Landau (<fr 2 ) 2 field theory in the outer space of the particle with the order 
parameter field <fr vanishing on the particle surface (see Ref. 1 19| , I, and Sec [II A[) . This 



allows us to resort to basic field-theoretical tools such as the renormalization group and 
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short- distance expansions which turn out to be extremely useful for the understanding of 
the polymer conformations in the presence of the particle (s). 

In the following we summarize our main results starting with the case of a single 
particle. The evaluation in I of the solvation free energy for immersing the particle in a 
theta solvent (i.e., ideal chains) has been generalized to the generic case of a good solvent 
(i.e., chains with EV interaction) by calculating the universal scaling function Y^d (1Z x /R) 



(see Eq. ( |1.5|) ). For estimates based on a systematic perturbative approach it is very 
useful to introduce the particle shape of a 'generalized cylinder' (see Eq. (|1.3|) ) which 
is characterized by the space dimension D and an internal dimension d encompassing 
cylinder, sphere and wall as special cases. The general results for Y^d(x) to first order in 
e = 4 - D are given by Eqs. (gjj) and (ggg ). 



(1) Our investigations in Sec. [IIB| of generalized cylinders with small curvature, i.e., 
R ^> 1Z X , provide strong evidence for the validity of the local and analytic Helfrich- 
type expansion conjectured in Eq. (|1.7|) . With the help of Eq. ( |2.39|) this expansion can 
be generalized to arbitrary spatial dimensions D so that we were able to obtain explicit 
expressions for the universal coefficients Act, A«i, Ak 2 , and Akq appearing in the Helfrich 
Hamiltonian to first order in e = 4 — D. While the results for the spontaneous curvature 
energy Ani in Eq. ( |2.42|) and the mean and Gaussian bending rigidities Ak 2 and Akq 
in Eqs. ( |2.43| ) and ( |2.44| ) are new, the result in Eq. ( |2.41| ) for the surface tension Aa has 
implicitly been noted before (see Eq. (4.7) in Ref. p3fl ). All coefficients have absolute 
values smaller than those of their ideal chain counterparts. The latter are given by the 
above expressions for e = 0. The decrease of the depletion effects due to the EV interaction 
can be traced back to a corresponding behavior of the profile Me of the end density (see 
Eqs. ( |2.1| ) and ( |2.2| )). The simplest case is the surface tension Aa which follows from the 
profile Me near a planar wall and for which the decrease is consistent with a corresponding 
decrease [19(a)] of the surface exponent oe in the behavior Me ~ (z/1Z x ) aE for distances 
z from the wall much smaller than TZ X . 

(2) For small particle radius, i.e., R -C TZ X , our results for Y^d i x ) t° first order in 
e confirm the validity of the power law ( |1.8[ ) within the region ( |1.9j ) and allow us to 
determine the e-expansions of the universal amplitude Aj,d (see Eq. fl2.52j )). The region 
( |1.9|) is shown shaded in Fig.|2| and includes the interior point {d,D) = (2,3) which 



represents a cylinder in three dimensions. This is different from the case of ideal chains in 
which Eqs. ( |1.8| ) and ( |1.9| ) are not valid below and on the line d = 2. Reliable estimates 
for the amplitudes and A 2 ^ corresponding to a sphere and a cylinder, respectively, 
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for chains with EV interaction in three dimensions have been obtained from the plausible 
assumption that the amplitude A^d as function of d and D forms a regular surface over 
the base plane (d, D) (see Fig.|||). The combination of the value of ^2,2 corresponding to 
a disc in two dimensions (see Table |) with the e-expansions of A^d i n Eq. (|2.52|) leads to 
the estimates in Eqs. ( |2.54| ) and ( |2.55| ) for A 33 and A 2 , 3 . 

(3) Estimates of the full scaling functions Y 3>3 (x) and Y 2 ^(x) for the solvation free 
energy of a sphere and a cylinder in three dimensions are shown in Fig. ^ in terms of the 
functions 63,3(0;) and 02,3(2) defined in Eq. (|2.56|) . This shows the crossover from the 
small curvature regime 1 < 1 with the coefficients 0(0), 0'(O), and 0"(O) of their Taylor 
expansions about the regular point x = being simply related to the surface tension 
Act, the energy Ak% of spontaneous curvature, and the bending rigidities Ak 2 and Akq, 
respectively (see Eqs. (|2.40f) and ( [2.561 )), to the small radius regime x ^> 1 with the 
power law 0,2,3(0; — > 00) — > Ad^x l ' v ~ x . As expected the curves in Fig. (| for chains with 
EV interaction are below the corresponding curves for ideal chains and imply a smaller 
solvation energy. For chains with EV interaction the exponent \jv — 1 is not a positive 
integer and the expansion of ©^.3 (x) or Y^ 3 (x) about x = cannot be a polynomial with 
a finite number of terms. This is in contrast with the solvation free energy of a sphere in 
a solution of ideal chains in which case Q(x) is a linear function of x (see Ref. Jll| or I). 

We continue by summarizing our results for the interaction between particles with 
a small radius. Since the values in Eqs. ( |2.54j ) and Q2.55p of the amplitudes A 3)3 and 
A 2)3 completely determine the Boltzmann weight in Eq. (|1 . 1 1|) of a small sphere and a 
thin cylinder, the interactions of these particles with other distant particles or walls are 
completely determined, too p6| , |38| , p9| , |57| , |64]| . 

(4) We have studied the interaction between a wall and a long thin cylindrical particle a 
distance T> apart with radius R and length I for the case R <C 7Z X , ^ <C I (compare Fig. |5|). 
The dependence on D of the polymer mediated free energy of interaction is proportional to 
that of the monomer density of a dilute solution of chains in the half-space without 
the particle (see Eq. ( |3.1| )). The same applies for a small sphere near a wall (compare I). 
Since has a point of inflection at T> ~ 1Z X the attractive mean force between a wall 
and a thin cylinder or between a wall and a small sphere somewhat surprisingly passes 
through a maximum as T> increases. The increase ~ X> l / U ~ l of the force per unit length 
l D ~ d and unit bulk pressure n p k B T with the distance V in the region i^cP C TZ X is a 
consequence of its length dimension d — 1, its independence of TZ X , and of the fact that 
the particle radius R enters the force only in the form of the power law R d ~ x l v according 
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to Eq. ( |1 . 1 ip . Our study of the situation of long chains is complementary to that of short 
chains, i.e., 1Z X <C R, considered in Ref. J5U| . In the latter case the attractive mean force 
of depletion is monotonically decreasing as T> increases. 

(5) The interaction between two small spherical particles A, B of equal size and with 
a distance tab = T> + 2R between their centers has been studied in Sec. [IIIB| both for 
R <C T>,7Z X and for R, V <C 1Z X . In the former case we use Eq. ( |1.16a| ) expressing the 
interaction in terms of ^4_d,_d, R, and the universal monomer density correlation function 
C2 of a single chain in unbounded space. In the latter case the 'dumb-bell' composed of the 
two spheres is small on the scale of TZ X and can in leading order be considered as a pointlike 
object. This gives rise to an expansion similar to Eq. (|1 . 1 1| ) in conjunction with the lower 
part of Eq. ( |1.12p in which, however, the amplitude Ad,d is replaced by an amplitude 
function AdAp^ depending on R and T>. This is another type of a short-distance like 
operator expansion which can be used not only for the effective free energy of interaction 
but also for other polymer properties — such as the monomer density profile — induced by 
the two spheres. Both cases overlap in the region fl«P< 1Z X in which the interaction 
free energy f AB per unit bulk pressure n p ksT is given by Eq. ( J3.14 ) . Numerical values 
of Ad jy and the universal bulk amplitude a in Eq. Q3.14|) are summarized in Table | for 
various space dimensions and both for ideal chains and chains with EV interaction. The 
value for a in D = 2 derived in Appendix is a new result for a self-avoiding chain 
in the unbounded plane. For D = 3 and ideal chains we explicitly calculate the two 
functions Ci and A4Ad,d (see Eqs. ( |3.10|) and (|3.20|) ) and thus present a complete and 
explicit expression for the free energy of interaction between two small spherical particles 
to leading order in the small quantity R/1Z X . In contrast to the polymer mediated force 
between a small sphere and a wall, for two spheres of equal size the force is monotonically 
decreasing in the whole range of T>. For the case of two touching spheres and arbitrary 
values of R/TZ X we consider an approximative form of f^ B (compare Eq. (|3.22|) ) and 



compare it with the results |L2] of simulations. 

(6) As an illustration for the nonpairwise character of the depletion interaction between 

particles we have evaluated an explicit analytic expression for the three particle interaction 

(3) 

JXbc m the cas e of small spherical particles and ideal chains in three dimensions. The 
expression follows by inserting the triple correlation function in Eq. (|3.15 ) of the monomer 
density in the unbounded solution in Eq. (|1.16b| ) and using that in this case D — 1/v — 1 
and Ao,d = A3 = 2tc. The result is valid in the region R rij,TZ x with rjj denoting 
the relative distances tab, fAc, or tbc between the spheres and is complementary to 
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the three-body results presented in Ref. |12| with R of the order of 1Z X . In order to 
convey an idea of the relative importance of one-, two-, and three-particle contributions 
we summarize the results 



(' 



{1 \ fl%, fi%,c) (4-1) 



= TkR-R* f 1, - 2 JL , 2R> r ^E±l^±l£A) 
V tab r AB r B crcA ' 

for the special case R <C <C 1Z X . For three small spheres configurated on an equilateral 

(3) (2) 

triangle with edge length r the interaction /X bc * s related to fX b f° r two spheres at a 
distance 2r via 

(Vile) /(-fi%) , = ( 4 - 2 ) 

This relation holds for an arbitrary ratio r/7Z x provided R Cr, 7^. 

Another interesting type of three-body depletion interaction arises for two spherical 
particles near a planar wall. If their radii are small this situation can again be system- 
atically investigated by means of Eq. fll.ll ) and the lower part of Eq. ( |1 . 1 2| ) . In the same 
spirit the investigations of three-body interactions could be supplemented to cover cases 
in which the distance between two of the spheres (or between one of the spheres and the 
wall) becomes of the order of R or smaller by means of the 'small dumb-bell' expansion 
(or an expansion which applies to a sphere close to a planar wall |p?|). 

Finally, we summarize some of the field-theoretic developments on which our treatment 
of the particle-chain interaction is based. 

(7) After a brief outline of the polymer magnet analogy in Sec. |II A| we relate the 
density profile Me of chain ends to the local susceptibility in the corresponding magnetic 
system for a generalized cylinder K in a dilute polymer solution (see Eq. (|2.27|) ). For such 
nonadsorbing chains the corresponding order parameter field $ vanishes at the surface 
of the particle. With the Gaussian order parameter correlation function outside K as 
the unperturbed propagator we use renormalized perturbation theory with respect to a 
(<fr 2 ) 2 interaction in order to obtain a systematic expansion in the EV interaction of the 
polymer quantities below the upper critical dimension D uc = 4. The behavior of our one 
loop expressions (see the function Cdij) in Eq. ( |2.31| )) in the limits corresponding to large 
R and small R is discussed in Appendix 0. 

(8) We verify to first order in the EV interaction that the same small radius amplitude 
appears for different properties of a generalized cylinder with a small radius R. In Ap- 
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pendix[FJwe write Eq. Ql.llj) in terms of fluctuating densities (operators) in the equivalent 
field theory. The universal small radius amplitude A^d f° r polymers is obtained from a 
corresponding critical amplitude A^d in the field theory by multiplying with a univer- 
sal noncritical bulk amplitude. In the two-point correlation function with distances of 
the two points from the generalized cylinder much larger than R there appears the same 
amplitude A^d a ^ the critical point of the field theory — where the correlation length 
£ + is infinitely large — as in the behavior of the field-theoretic excess susceptibility of 
the generalized cylinder for £+/R ^> 1. The latter is related to the power law behavior 
flL8|) of the function Y dD (x) for x = 7Z X /R ^> 1. These considerations are important 
to understand that the mechanism behind the small radius expansion is basically of the 
same type as that behind the well-known short distance expansions in field theories with- 
out boundaries f39,64j|. Moreover, in case of a sphere our result for Ad,d to first order 
in e confirms that this amplitude can be reduced to bulk and half-space amplitudes as 
predicted from a conformal mapping |65| (see the last but one paragraph in Appendix[B|). 

(9) By studying the energy density profile in a Gaussian field theory with boundaries 
we explicitly verify that not only a single sphere but also a 'dumb-bell' composed of two 
spheres of equal size can be considered as a pointlike perturbation on sufficiently large 
length scales. At bulk criticality the profile for the dumb-bell can be obtained by means 
of a conformal transformation from the known profile between two concentric spheres. 
For ideal polymer chains we thus find the explicit form (see Eq. ( |3.20| )) of the amplitude 
function A4Ad,d addressed in paragraph (5) of this Summary. 



ACKNOWLEDGEMENTS 

We thank T. W. Burkhardt for helpful discussions. The work of A. H. and S. D. has 
been supported by the German Science Foundation through Sonderforschungsbereich 237 
Unordnung und grofie Fluktuationen. 



48 



APPENDIX A: THE FUNCTION C d {r) 



The results of Sec.[TJ are based on the behavior of the function Qd,D (v) m Eq. ( |2.30|) , in 
particular on the behavior of Cd(rj) in Eq. (|2.31|) . The difficult part of the corresponding 
calculation consists in performing the sum over n and the double integral over q and ip 
in order to calculate C<f(V) according to Eqs. ( |2.15b| ) and ( |2.31b| ). Here we derive the 
asymptotic expansions of Cd(r) for large and small r, respectively, and give numerical 
values of Cdij) for the crossover region < r < 3. 



1. Cd(j) for r — > oo 



We calculate the coefficients C[ and in Eq. fl2.32| ) for d — D, 3, and 2 by expand- 



ing the rhs of Eq. (|2.31b| ) for large r. To this end we need the behavior of the integrand 
in Eq. ( [2.31b|) for R/jy/t = \fr large and (y± — R) [i\ft — — 1)^/t = s arbitrary. This 
is consistent with the expectation that for the small curvature expansion the important 
regime in terms of polymer variables is R/1Z X large and (y± — R)/1Z X arbitrary. 

(a) d = D: Since and C 2 belong to the one loop contribution of Qd,d{v) we 
need to consider only C/ and C^ 1 (compare the remarks below Eqs. ( |2.20| ) and ( p.31b|) ). 
The central part of the calculation consists in expanding g s (ip,r,e = 0) in the integrand 
on the rhs of Eq. ( |2.31b| ) for r large and s arbitrary. Since Wn^(0) = (n + l) 2 /(27r 2 ) for 



a = 1 in Eq. (|2.15b| ) the quantity g s (ip, r, 0) is, apart from a factor — ip 2 / (2n 2 ) , given by 

±« 2 w£r) [K " (3+v ~ T)? - (A1) 



n=0 



A first hint on how to evaluate the sum ( |A1| ) for large r can be gained from recognizing 
that its leading behavior corresponding to a vanishing curvature must describe the half- 
space bounded by a planar wall. This is discussed after Eq. (|2.15c|) and shows that the 



ratio n/R has the meaning of the length of a wavevector parallel to the wall and that all 
values of n are important for which [n / R) / (fiVt) = n/ ^fr = wor [n/ R)(y± — R) = ns/ 'y/r 
are of order unity. Thus for the general expansion for large r a large number of terms will 
contribute and the sum can be replaced by an integral plus corrections according to the 



Euler-MacLaurin formula 66 



oo 00 

£ F(n) = [ dnF(n) +±F(0) - |f F'(0) - |i F"'(0) + ... . (A2) 
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Here are Bernoulli numbers and the function F(n) can be read off from the expression 
(|AT|). For case (a) the analysis of this expression shows that all contributions on the rhs 
of Eq. (|A2j ) apart from the integral lead to orders of r -1 / 2 higher than needed for the first 
three terms on the rhs of the expansion (|2.32|) of C^(t) (but compare case (b) below). 
Upon introducing uj instead of n as the integration variable the expression ( pTT| ) turns 
into 

oo 

j3l2 r^Mt^)]' (A3) 

J K a {a/uj) L V uj^/t /J 

o 

with a = Uy/r. For large r the integral ( |A3| ) can be simplified by employing the uniform 
asymptotic expansion for large orders a of the modified Bessel functions I a and K a which 
is provided, e.g., in the sections 9.7.7 and 9.7.8 of Ref. [40(a)]. In addition to the leading 
term (compare the discussion after Eq. (|2.15c|) ) now also the correction terms containing 



the functions Uq, ui, and Ui given in section 9.3.9 of the above reference have to be 
included. By inserting this simplified integral into g s in Eq. ( |2.31b| ) one finds that the 
first three coefficients on the rhs of the expansion (|2.32|) of C^{r) are determined by a 



number of double integrals over s and uj which can all be calculated in closed form. This 
reproduces the expression ( |2.35|) for C — and thus checks the assumption leading to it — 
and yields the expressions for and C 2 in Eq. ( |2.36| ). 

(b) d — 3: Due to the additional integration over q in Eq. ( |2.15b| ) the expression 
corresponding to ( |AT| ) now reads Y2^=o F(n+l/2) where, using wi 1/2) (0) = (n+l/2)/(27r) 
and substituting k = qr~ x l 2 in Eq. ( |2.15b|) , 



OO 



Fin) = rm\ ^|r^=r \KJ(s + V^)^+l)} 2 . (A4) 

J K n (\/T\/K, 2 + l) L \ I- 







From the Euler-MacLaurin formula ( |A2| ) one infers that, in contrast to ), apart 

from the integral on the rhs also the terms proportional to F(l/2) and to F'(l/2) have 
to be included in order to obtain the first three terms on the rhs of the expansion (|2.32|) 
of Cs{t). Proceeding in the same way as in case (a) one is led to consider modified Bessel 
functions I a and K a with order a = ujy/r\/ k 2 + 1 and triple integrals over s, u, and k. One 
reproduces again the expression ( |2.35| ) for Co and finds, using B 2 = 1/6, the expressions 



for C/ 3) and C 2 (3) in Eq. (fZ37|). 

(c) d = 2: In this case the procedure is quite similar as in case (b). The expression 
corresponding to ( |A1| ) now reads F(0)/2 + Yl^Li F(n) where 
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F(n) 



dn K 



K„ 



tVk 2 + 1 



(A5) 



The analysis shows that only the integral on the rhs of the Euler-MacLaurin formula ( |A2| ) 
contributes to the first three terms on the rhs of the expansion fl2.32p of ^(r) (compare 
cases (a) and (b) above). One finds again the expression fl2.35|) for C and in addition the 
expressions for C} 2) and C 2 (2) in Eq. flpl). 



2. C d (r) for r ^ 



The leading behavior of C d (r] — > oo) in Eq. ( [2.471 ) can be inferred from the behavior 
for r — > of the quantity 

W = ~f^ C d (r) (A6) 

with C d (r) from Eq. ( |2.3 1| ) . The behavior of 1 d {r — > 0) exhibits two types of leading 
terms. The first is the logarithmically divergent contribution — a/(47r 2 ) ln(l/-^/r) which 
follows from the behavior g s (ip, r, 0) — > — a/(47r 2 ) for 1 <C t/> <C 1/ \fr as mentioned below 
Eq. (|2.15d|) . The second contribution is independent of r and requires special care. Its 
evaluation is facilitated by splitting T d {r) according to 

Xdir) = H d (r) + J d (r) , (A7a) 



where 



7~Ln 



T 



dijjijj' 



(A7b) 



x 



1- ip- 



K a (^) 
K a (^F) 



g s (iP,r,0) - sj as) (^,0) 



Mr) 



(A7c) 



Here we have used Eqs. ( |2.9| ) and ( |2.12| ) and we have added and subtracted the function 
gs (ip\/T, 0) which is defined as in Eq. ( ^.15d| ) and represents the behavior of g s (ip, r, 0) 
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for 1 < if), ^~ 1/2 - In H d one can interchange the order of the integration over ijj and the 



limit r — > |67| which results in the finite limit 



oo 

?k(T-0)-B d = I [1-^ 2 1 2 7,(^,^ = 0) + 



(A? 



where the function 



lMi £ ) = 9s(ip,T= 0,5) 



(A9) 



can be read off from Eq. ( |2.15b| ). The integral in Eq. ( |A8| ) is well-defined since 7 s (?/>,0) 
tends to — a/(47r 2 ) for large ip so that the logarithmic singularity is removed. The integral 
in Eq. (|A7cj ) can be carried out explicitly and leads in conjunction with Eqs. ( |A"7| ) and 
to 



1 d {r -> 0) -> B d + 



47T 2 



In 



r 



#(d/2) C E 
ln2 + 1 -^ + f 



(A10) 



where \l/ is the psi-function and Ce denotes Euler's constant. Inserting Eq. (|A10|) in 
Eq. ( |A6| ) and carrying out the inverse Laplace transform in Eq. ( [2.31a| ) leads to the result 
for Y d ^(x — > 00) in Eq. (|2.48|). We conclude this subsection by calculating the number 



B d for d = D, d = 3, and d \ 2 (see Eq. flTSOD). 

(a) d = D: Since Xd belongs to the one loop contribution of Yd,d we need to consider 
X4 only (compare the remarks below Eqs. ( |2.20| ) and ( |2.31 bj) ) . This amounts to inserting 
a = 1 into Eq. ( |A5| ) and the function 7 S corresponding to a sphere in D = 4 which is 
given by 7,(^,0) = -(4tt 2 )- 1 [1 - ^~ 2 ]~ 2 |6| yielding i3 4 = 0. 

(b) d = 3: For 2 < d < D the quantity B d does not vanish and can be evaluated 
numerically. For d = 3 this leads to the value for B% given in Eq. (|2.50 ). 

(c) d \ 2: In this limit £>d can be calculated exactly. It is useful to substitute a = ip 2a 
in Eq. (|A8|) and to carry out the limit a — (d — 2)/2 \ for fixed a in the ensuing 
integrand. One finds that only the term for n = in Eq. ( |2.15b| ) survives this limit with 
the result 



B 2 



da a 



-1 



a 



a 



8tt 2 



(All) 
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3. Cd{r) in the crossover region < r < 3 

For the convenience of the reader in Table [TI| we give some numerical values of C^{r). 
From these values an approximation for the full function Cd(r) can be constructed by 
using its asymptotic behaviors for r —>■ oo and r — ► as derived in the above subsections 
and by appropriate interpolation. 



TABLE II. Numerical values of In C d {r) (see Eq. ( ^31b| )). 



lnr 


d = 2 


d = 3 


d = A 


-10 


18.816 


21.308 


22.223 


-9 


16.795 


19.175 


20.108 


-8 


14.773 


17.027 


17.980 


-7 


12.755 


14.865 


15.835 


-6 


10.744 


12.692 


13.672 


-5 


8.748 


10.511 


11.488 


-4 


6.773 


8.328 


9.282 


-3 


4.830 


6.160 


7.057 


-2 


2.931 


4.024 


4.836 


-1 


1.085 


1.946 


2.640 





-0.700 


-0.054 


0.504 


1 


-2.422 


-1.966 


-1.540 



APPENDIX B: SMALL RADIUS EXPANSION TO ONE LOOP ORDER 

The relation ( |1 . 1 1| ) for polymers is — via the polymer magnet analogy — closely related 
to a corresponding small radius expansion (SRE) in a ($ 2 ) 2 field theory with the Boltz- 
mann weight exp(— A7ix{^}) which describes the presence of the generalized cylinder K 
(compare Sec. II A| and Appendix 0). Here we shall illustrate the SRE by considering the 



two-point correlation function at the critical point of the field theory in one loop order. 
This is particularly well suited to reveal the mechanism behind the SRE. Moreover it 
provides a significant check for the operator character of the expansion because we shall 
find the same small radius amplitude A djD as in Sec. |II Q . 
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Keeping u and e = 4 — D as independent variables the SRE can be written in the form 
exp(-AH K ) oc 1 - T(nR, u; e, d) /i 2 ~ d Z t u K + ... (Bl) 



with 



/ d s r\\ $ 2 (r ± = 0,r||) , d < D 



$ 2 (0) 



d = D 



(B2) 



Here /i Z t ujk is a renormalized and dimensionless operator and 

T{^iR,u;e,d) = — Ajp (fiR) d ~ 2 [1 + uFi(nR; e, d) + 0(u 2 )} 



(B3) 



has an expansion in terms of u with the coefficient — A$ = 2ir d / 2 /T(a) = aVt d of the 
leading term corresponding to the Gaussian model (see Eq. (4.6) in I). The functions Fj 
can be expanded in terms of e with coefficients which depend on nR only via powers of 
In(fiR). In particular we shall find from the critical two-point function that 



Fi{jj,R\ e, d) 



H+2 



4tt 2 

In(fiR) +fi + e d + B d + 0{e) 

a 



where 



ed 



1 + 



lnvr *(d/2) 



(B4) 



(B5) 



the quantity B d has been introduced in Eq. (|A8[). The ellipses in Eq. ( ]B~1"| ) stand for 
contributions in which higher powers of R are multiplied by powers of lni?. Standard 
renormalization group arguments imply that for large fiR the function T is proportional 
to R d ~ l/u and that the rhs of Eq. can be written as 1 + ArR 1 ' v ^>k where 



A K = fi 2 - 1/u Z t D L (u)F(l,u*;e,d) 



(B6) 



with Dl from Eq. (|2.28 ). The universal polymer amplitude A^d in Eq. ( |1.11| ) is related 
to A K = A K {N) via m 



A 



d,D 



A K (0)2 I ,- 2 Z- X LTZ^ . 



(B7) 



By using Eq. fl2.28|) one finds that the nonuniversal quantities /x, Z t , D L , fi cancel and 

-4vr 2 3 In 2 *(d/2) 



2lx d/2 



r(a) 



4la Bd+ 2 2 



0(e 2 )} 



B8 
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which indeed reproduces the first-order e results of A^d in Eq. Q2.52 ). 



We now verify Eqs. ( |B1[ ) - ( |B5[ ). Consider the two-point correlation function 



($j(r)$fc(r ')) of the field theory described by Eq. (|2.6|) at its critical point. For u = the 
SRE follows from the explicit expressions in Eq. (|2.7Q for the Gaussian propagator which 
by using Wick's theorem lead to 

(a*$ i (r)$ fe (r')) M o] (B9) 
= (AP)- 1 limi?- 2a {($ J (r)<|. fc (r')) [ o ] -(^(r)$ fc (r'))M 

-ft — >U K 



Here ( ) is a cumulant average with the subscript [0] indicating u = and with b denoting 
the unbounded bulk space in absence of K. Obviously Eq. (|B9|) verifies the SRE for the 
Gaussian model. 



Consider now the first order in u contribution: 

($,(r)$ fe (r')>[i] = -5 lk ^^8n 2 f^uR 2a J(v,r') (BIO) 

where 



J(r,r') = J d D yG(r,y;R)G(r', T ,R)I(y x ,R), (Bll) 



y±>R 



I(y ± ,R) = R- 2a G(y,y;R) = — r { ($ 2 (y)) [0] - ($ 2 (y)) M o]} • (B12) 



The first order expression given in Eq. ( |B10| ) has the same structure as the one in Eq. ( |2.11| ) 



and we have used Eq. (2.16a) . Note that in the present dimensional regularization scheme 



and at t = the bulk quantity G&(y, y) = ($ 2 (y))fe ) [o]/W vanishes. We have exploited 
this in order to write the last expression in Eq. ( P12|) in such a form which allows us to 
make contact with Eq. ( |B9| ) and which implies 



l(y±,0) = ^(^# 2 (y)) M0 ] = -yl d+£ ^{l + ee d + 0(e 2 )} . (B13) 
The function I(y±, R) is related to 7 S (^,£) in Eq. ( |A9| ) by 

yt £ l(y±,R) = ls(y±/R,e) (B14) 
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and Eq. ( |B13| ) is consistent with 7 s (oo,0) = — a/(4iT 2 ) as mentioned below Eq. (|^9|). In 
order to verify Eqs. ( B~T| ) - flB"5] ) we decompose J(r, r ') according to 



J — J(i) + J(u) + J(in) 



(B15) 



with 



^)(r,r') = J d D y G b (r, y) G b (r ', y) I(y ± , R = 0) , 



(B16a) 



J ( «) (r, r ') = - y tPy G 6 (r, y) G b (r y) I(y u R = 0) , (B16b) 

0<y x <-R 



■/(«*) (r, O = J d D y {G(r, T ,R)G{r',y;R)I(y ± ,R) 

-G b (r,y)G b (r',y)I(y L ,R = 0)} . 



y±>R 



(B16c) 



In the following we analyze the behavior of the rhs of Eq. ( B10|) for R — * which arises 
from each of these contributions J^, J^), and 

Rewriting /(t/j_,0) in the integrand on the rhs of Eq. ( |B16a| ) by means of the first 
equality in Eq. ( pi3|) one finds (see Fig.||) 



(^■(r^rOWw - APR 2a (u K ^(r)Mr')) b , m ■ 



(B17) 



In order to evaluate the leading contribution of Jiu) for small R one can set y^ = in 
the two bulk propagators G&(r, y) and G b (r ', y) in the integrand on the rhs of Eq. (|B16b| ). 
Its remaining dependence on y± as given by the last expression in Eq. ( |B13| ) leads to a 
pole in e. This results in (see Fig. |]) 



(^•(r)$ fe (r ')}[!], («) - 4 0) iJ 2a K$ 3 (r)^(r')), 
x {(Z t )[i] + u \n{nR) + fi + e d } 



(B18) 



where {Z t )[x\ = (N +2)u/ (3e) is the contribution to Z t of first order in u (see Eq. ( p,16b|) ). 
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• • • 



CO 



X = l 




K 



Eqs. (B.18), (B.22) 



Eq. (B.17) 



FIG. 6. Diagrammatic representation of (lok ^i( r )^i( r '))b, [A] appearing on the rhs of 
Eqs. (B17), (B18), and ( B22| ). The solid lines correspond to the bulk Gaussian propagator 
[o] = Gb and A = 0, 1 denotes the loop order. The wiggly lines indicate the insertion of the 
operator u>k located at the 'axis' of the generalized cylinder K. 



For the leading contribution of J(iu) for R — > it is sufficient to confine the integration 



over y± on the rhs of Eq. ( |B16c|) to the restricted region R < y± < \l Rr± with r^ <] 



min (r±,r±). The reason is that in the remaining integration region not only the ratios 
R/r± and R/r'^ but also R/y± are small so that in leading order one can insert R — into 
the first term in curly brackets, which is then cancelled by the second term. In particular, 
in the restricted region y± is smaller than r± and r ± . By inserting the representation 
( p.7| ) for the external legs G and Gb in Eq. ( |B16c| ) one finds that only the terms for n = 
contribute to the leading behavior of Jim) for which one obtains 

J(^)(r,r') -> ^R e (aj K ^ 1 (r)^ 1 (r , )) bt[0] (3 d , (B19) 

oo 

(3 d = J #^" 1+£ { [l-^- 2a ] 2 -f s (^,e) -7.(°o,e)} (B20) 
i 

with 7 S from Eq. (|B14j ). The quantity (3d arises as the limit for R — > of the expression 
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x 



d^iP~ l+e [T{a+ l)] 2 



-2a 



\ 2 
I a {RP) 



(B21) 



K a (Rpy 



K a (RP^ ls (^,e) - [I a (RPiP)] 2 ls (oc,e)} . 



Of course, [3d — > Bd for e — > and the present procedure for R/r 



(<) 



at the critical 



point of the field theory leading to Eq. (|B20|) should be compared with the procedure for 
R/£+ -> leading to B d in Eq. (TAE|) . Equation ([BT9[) implies (see Fig.|) 



($,(r)$ fc (r')) W 



(iii) 



„4° } i? 2 ° (^$,(r)$ fc (r')) M 



TV- 



(B22) 



2 4vr 2 

-u fi d + 0{ue) . 

a 



Equations (|Bl7| ), (gig ), and ( p22|) corroborate the SRE in Eqs. flBl]) - flB5|) to first order 
in m in the case of the critical correlation function. Note the recurrent character of the 
SRE which is typical for operator product expansions [39(b)]. A graphical representation 
of the bulk correlation function with insertion of the operator ujk is shown in Fig. 

Apart from particles in a polymer solution there are other physical systems the SRE 
can be applied to such as spherical or cylindrical particles in liquid 4 He near the A point 
and nonmagnetic inclusions in a ferromagnet of Ising or Heisenberg type near the Curie 
point. For these systems the parameter Af takes the values 2, 1, and 3, respectively. A 
useful characterization of the small sphere or the thin cylinder in these cases is provided 
by the universal amplitude 



A d , D {U) = -A K {M) 



with the amplitude B<$>2 of the bulk correlation function (<fr 2 (r)<fr 2 (0))fc = B^r 



(B23) 



at the critical point. For example the change in free energy per unit 'length' I s which 
arises upon immersing the generalized cylinder K into the bulk system displays a singular 
dependence on t ~ (T — T c )/T c given by 

k B T 



(B24) 



sing 



-k B T c R d -^A K (^) [(3>\tUn 9 
k B T c R d ^» C {D - 1M A d , D {M) E(N) 



with the universal bulk amplitude 
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E{N) = l($XtUn g t; D - 1/v (B&/M)-V 2 (B25) 

which characterizes the temperature dependence of the bulk energy density. From 
Eqs. (|B~1| ) - ( |B6| ) and from the dependence of B§2 (A/") on e one obtains the following 
explicit expressions: 

A D)D {M) = + 0(e 2 ) , (B26a) 
V 2 

1 3 ,d(AT) = JL (i + I [c B + hi7r (B26b) 

a/2 7r *• 2 L 

JV+2 



+ ^rl( 167r ^ 3+21n2 - 1 )]} + G(£2) ' 



and 



A 2 , D (AT) = e^^2 3 /\l3 2 + 0(e 2 ) . (B26c) 



The first-order e result (|B26a| ), which we have obtained by carrying out the calculation 
directly in the outer space of a sphere, confirms the prediction 



A D , D (M) = yJ(A$ 2 ) 2 /(AfB^) (B27) 

which follows from relating the half-space (hs) profile (® 2 (z)) hs = A$\2z)-( D -^ with 
the distance z from a planar wall with Dirichlet boundary conditions O at the bulk critical 
point to the profile (<fr 2 (r)) in the outer space of the sphere by means of a conformal trans- 
formation |55| (for M = 1 compare the explicit result in the first Eq. (20) of Ref. [36(a)]). 
The consistency of the above results is expected but remarkable since the finite conformal 
map changes the geometry under consideration. 

It is helpful to summarize the relationship between the universal amplitude A^d for 
polymers and the universal amplitude A^d for the field theory in terms of the symbolic 
equation 

A diD V = A d>D y/N/Bv $ 2 (B28) 

which applies inside averages or correlation functions for M \ with \1/ defined in 
Eq.(p7|). Since * = y/b*/B& $ 2 with = {Jll /v /{2L )fB^2 from Eq. flCIDD one 
has 
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A d , D = (A d , D ^N/b* ) . (B29) 

For a spherical particle, in particular, the polymer amplitude Ad d can be expressed in 
terms of the critical universal amplitude ratio on the rhs of Eq. ( |B27|) and the noncritical 
universal bulk amplitude hq,. In D = 2 both are explicitly known [19(b)] and lead to the 
value A 2> 2 = 3.81 in Table |. 

APPENDIX C: SHORT DISTANCE AMPLITUDE FOR POLYMER DENSITY 

CORRELATIONS 

Here we calculate the universal amplitude a in Eqs. ( |3.7| ) and ( |3.8| ). While for D = 4 
it coincides with the corresponding ideal chain value 

a (id) = n -2 from Eq m D = 3, 

2, and 1 the amplitude a is different from a^ d \ 

(a) D = 3: In this case results are available |T8| , |53|| for the normalized scattering form 
factor H(Q) which is defined for general D by 

C 2 (r,0) = TZ^ J 0^ exp(*Qr)iJ(Q) , (CI) 

R D 

where H(0) = 1 as implied by Eq. ( |3.5| ). The amplitude hoo in the power law 

H(Q ^ oo) ^ hoc {Q 2 K 2 J2)- 1 I^ (C2) 

is related to a by 

a = hoo k-W r (^^) jv (i) . (C3) 



From the accepted |T8| , |53l approximate value /loo « 1.1 in D = 3 one infers via Eq. ( |C3| 



the value a w 0.13 (see Table |). 

(b) D = 2: In this case one can obtain a fairly accurate estimate for a by combining 
a numerical estimate for a ratio of gyration radii of ring- and open-chain-polymers with 
conformal invariance and Bethe ansatz results for the 0(N) vector model by invoking the 
polymer magnet analogy (PMA) |15|,i~7||. By using the language of the Ginzburg-Landau 



field theory (compare Sec. [11 A| ) the necessary relations of the PMA can be written in a 
way which makes the generalization to D = 2 obvious. The polymer average 
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{p{v A )p{v B )p{x c )} y (C4) 
_ £(tt(r A )tt(r B )tt(rc)$i(y) / d D y' ^(y')) 

£($i(y)/A'*i(y')) 

is expressed in terms of cumulant averages ( ) of the field theory. Here C = £(io — > Lq) 
denotes an inverse Laplace transform defined as in Eq. ( |2.5| ) and relates the strength to in 
the thermal perturbation 

H th = J d D rT(v) , (G5) 

T(r) = |d> 2 (r), (C6) 

of the Hamiltonian at the critical point of the field theory to the bare 'chain length' 
Lq which — apart from a nonuniversal proportionality factor — equals the number of 
monomers in the polymer chain. The scaling dimension of the quantity 

¥(r) =^^T(r) (C7) 

equals its inverse length dimension D — \ jv. Here E = t L is the exponent which 
appears in C (compare Eq. ( |2.5| )). The rhs of Eq. ( |C4| ) has the normalization property 
that by integrating the numerator over, say, ya one can replace J d D rA^(rA) by 
(compare the discussion related to Eq. (18) in Ref. |62]]). This is consistent with the 
corresponding normalization property J d D rA p{^a) = T^x for the lhs of Eq. ( |C4D as 
implied by Eq. ( |1.13| ). 



Short distance properties such as those in Eqs. (|3.7| ) and ( |3.8|) follow from the operator 
product expansion (OPE) 

*(r A ) *(r fl ) - ar,f lW * + r B )/2) , (C8) 

which is equivalent to the well-known OPE of energy density operators [39(b), 64]. The 
amplitude a is expressed as 



a = C V ■ (C9) 

Here bq, is the universal bulk amplitude in 

(*(r)*(0)) crit = bvr- 2 ^- 1 ^ (CIO) 
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with ( ) crit denoting the average at the critical point of the field theory and ( is the 
amplitude which replaces a in the corresponding OPE for the normalized energy density 
^ = \£?/^/%. By using results of Refs. 0101 one finds that in D = 2 for Af \ 

A/^C - (216 vr)^ ^ - 1.21 (Cll) 

(compare Eq. (7.161) in Ref. [19(a)] where ( is denoted by c 2 ). The amplitude bq, has been 
calculated in the Appendix of Ref. [19(b)] by using results of Ref. []73| so that in D = 2 
for Af \ one has 

/p. \ 2/3 

^- 1/2 Vh-^ 1 (\K/nL g ) (ci2) 



where k = 0.226630 and ftl/ft^ne ~ 6.85 |73j] is the ratio of 1Z% of an open polymer chain 



and the mean square radius of gyration 7£ r 2 ing = T^x, ring + ^y, ring °f a polymer with 



the same number of monomers. Equations (|C9|) - (p!2 ) lead to the value a ~ 0.278 in 
Table | 

(c) D — 1: In this case the behavior of a chain with excluded volume interaction is 
that of a rigid rod of length 1Z X . Thus v — 1 and C 2 (r, 0) equals — |r| for |r| < ^ 
while it vanishes for |r| > 1Z X . The assumption that Eq. (|3.7|) still holds for D = 1 leads 
to the value a = 1 in Table [D. 
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